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Abstract 



Quasi-one dimensional quantum spin systems have been one of the most active sub- 
jects of research in condensed matter physics in the last two decades. One of the 
main reasons is that such systems are simple model systems, which are mathemati- 
cally more tractable than higher dimensional systems, and at the same time exhibit 
a variety of novel magnetic phenomena due to strong quantum fluctuations. In this 
dissertation, experimental results on three quasi-one dimensional Heisenberg antifer- 
romagnets are reported. 

CuCb • 2DMSO (CDC) has previously been reported to be a spin-1/2 one di- 
mensional antiferromagnetic Heisenberg chain with J/ks ~ 17K in zero field. Field- 
dependent specific heat measurements have been used to explore the magnetic prop- 
erties of CDC. We have found an interesting phase diagram when the system is placed 
in an applied field. The 3D ordered phase is suppressed when the field is parallel to 
either the crystalline a axis or c axis, while the 3D ordered phase is enhanced when 
the field is parallel to the b axis. A field-induced gap is also observed when the ex- 
ternal field is parallel to the a and c axis and the magnitude of the gap changes with 
orientation and magnitude of field. The primary mechanism for the gap formation 
could be an effective staggered field due to both the alternating g-tensor and the 
Dzyaloshinskii-Moriya interaction. The magnetic structure in the ordered phase and 
the critical exponents of the 3D ordering phase transition have been determined on 
the basis of neutron diffraction measurements. 

Ni(C5Di4N2)2N 3 (PF 6 )(NDMAP) is a quasi-one-dimensional antiferromagnetic 
spin-1 chain. We have characterized the spin Hamiltonian of NDMAP by zero field 
inelastic neutron scattering. The intra-chain coupling constant, the single anisotropy 
parameters, and the Haldane gap values of the system are determined. In an applied 
magnetic field we found commensurate transverse magnetic ordering by means of 
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neutron diffraction. Depending on the direction of the applied field, the high field 
phase has either three dimensional long range order or quasi-two-dimensional short 
range order. We also performed inelastic neutron scattering experiments in a magnetic 
field below and above the critical field H c at which the gap closes to probe dynamics 
spin correlations. For H < H c , the behavior of NDMAP is well described by first 
order perturbation theory and is similar to that previously found in the extensively 
studied Haldane gap compound Ni(C 2 H 8 N 2 ) 2 N0 2 (C10 4 ) (NENP). For H > H c , we 
have observed quasi-elastic excitations with broadened q-dependence. Various models 
based on topological excitations with low energy phase fluctuations may account for 
the data, including a sine-Gordon model, have been shown to account for the phase 
diagram, and a free fermion model has explained the observed q-width. 

Y 2 BaNi05 is a quasi-one-dimensional spin-1 antiferromagnet, in which holes can 
be added by substituting Ca 2+ for off-chain Y 3+ . While the material remains insulat- 
ing, calcium doping introduces holes that reside on super-exchange mediating oxygen 
atoms. Thus, Y 2 _ a .Ca a; BaNi05 is a one-dimensional analog of the cuprates, where 
off-Cu0 2 -plane chemical impurities donate holes to the Cu0 2 planes. Inelastic neu- 
tron scattering experiments were performed to study the magnetic excitations in the 
pure, £=0.04 and x=0.1 samples. The parent compound is a spin liquid with macro- 
scopically coherent quantum ground state and no static order. The main feature of 
the excitation spectrum is an ~ 8 meV Haldane gap separating a triplet band from 
the singlet ground state. At energies above the spin gap, the triplet excitations of the 
parent compound persist with doping. Below the gap, new excitations appear with an 
incommensurate structure factor. The incommensurability saturates with increased 
doping. We show that the incommensurate peaks arise because of the characteristic 
spin density modulation that develops around a defect in the singlet ground state of 
a quantum spin liquid. The experiment provides a quantitative view of the magnetic 
polarization cloud associated with holes in a doped transition metal oxide. 

Work in this dissertation was performed under the supervision of the author's 
advisors: Professors Collin L. Broholm and Daniel H. Reich of the Department of 
Physics and Astronomy at Johns Hopkins University. 
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Chapter 1 
Introduction 

Although it is one of the oldest known and longest studied physical phenomena, 
magnetism continues to be an active and challenging subject. It is an experimentally 
driven field, yet it has also been a continuous source for new theoretical ideas. Much 
of current technology is based on magnetism, and future technical developments rely 
on further progress in understanding and controlling magnetism on the nano-scale 
ranging from spintronics to quantum computing. The interplay of magnetism with 
other areas of science and engineering has also stimulated a large number of research 
activities in recent years. 

One dimensional antiferromagnets have been the subject of intense theoretical and 
experimental investigations for several decades. One of the reasons is that they are 
simple model systems, which are well suited to test theories of many body physics. 
Furthermore, the effects of quantum fluctuations are more significant in one dimen- 



2 



sion, resulting in many interesting magnetic properties. On the other hand, recent 
developments in the synthesis of complex materials and in experimental techniques 
have reinvigorated the field from the experimental point of view. Moreover, many 
of the compounds that exhibit high-T c superconductivity are quasi-two-dimensional 
Heisenberg antiferromagnets. Antiferromagnetism in quantum spin chains has at- 
tracted great attention, as a possible one-dimensional analog to these systems. 

1.1 Conventional antiferromagnets 

The low energy properties of a non-metallic solid can be described by an effective 
spin Hamiltonian known as the Heisenberg model 



where J is the strength of the exchange interaction, S is the spin operator, and < ij > 
represents a pair of nearest neighbor (n.n.) spins. 

If the Si were classical vectors of length S, each n.n. spins pair has an interaction 
energy JS 2 cos#, where 9 is the angle between n.n. spins. If J < 0, the lowest energy 
is achieved when 9 = 0, i.e. the spins are parallel, giving rise to a ferromagnetic 
ordered ground state. If J > 0, the lowest energy is achieved when 9 = it, i.e. the 
n.n. spins are antiparallel; this ground state is called antiferromagnetic Neel order. 

Magnetism, however, is a purely quantum phenomenon. Spin operators satisfy 




(1.1) 



<ij> 
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commutations rules: 



[S a ,S p ] 




s 2 



5(5+1) 



(1.2) 



where S a are the vector components, and e Q/ g 7 is the Levi-Cevita tensor. Eq. 1.1 can 
be rewritten as 



where Sf = Sf ± iSf are the raising and lowering operators. 

For a ferromagnetic material, the classical ground state is still the quantum me- 
chanical ground state. However for an antiferromagnetic material, the Neel state is 
not even an eigenstate of the Heisenberg Hamiltonian because of the spin-flip terms 
Si + Sj~ and S^Sj + , and therefore can not possibly be the ground state. 

To study how much the Neel state deviates from the exact ground state of an 
antiferromagnet, a standard approach is linear spin- wave theory. In three dimen- 
sions, spin wave theory shows that the Neel state is a good approximation to the 
ground state of most three dimensional antiferromagnets. In one dimension, how- 
ever, the analysis reveals that S— < S z >, the deviation from maximum sublattice 
polarization, diverges due to quantum fluctuations, indicating that the Neel state is 
completely unstable no matter how large S is. In the following sections, we shall 
discuss the one dimensional spin system in greater detail. 




(1.3) 



<i,j> 



1.2 Half-integer spin Chains 

The 5 = 1/2 one dimensional spin chain is exactly solvable by Bethe Ansatz[l. 2]. 
A set of exact results, such as the eigenfunctions, and the ground state energy, have 
been obtained. The equal time correlation function is calculated as [3] 

(S -S 1 )^(-l)'|/|- 1 (ln/) 1/2 . (1.4) 

Because there is no characteristic length scale, this system is said to posses quasi- 
long-range antiferromagnetic order. 

The lowest-energy excited states have also been calculated by des Cloizeaux and 
Person (dCP) and the spectrum is given by [4] 

71" 

e d cp{q) = ^J\smq\, (1.5) 

where q is the wavevector transfer along the chain. The dCP result is actually the 
lower-bound of a continuum of excited states with an upper boundary of [5, 6] 

e u (q) = tyJ\ sin||. (1.6) 

It is illuminating to map the spin 1/2 chain onto the one dimensional system of 
interacting fermions[7]. For simplicity let us consider a one-dimensional XY spin 1/2 
chain in a magnetic field applied along z direction. The Hamiltionian can be written 

as 



(1.7) 



and it can be converted to the fermionic Hamiltonian by the Jordan- Wigner trans- 
formation [8, 9]. The transformation introduces a pair of fermionic operators a*, a/" 
on each site and expresses the spin operators in terms of fermions 



= af exp(i7r^a+o j )> (1-8) 



j<i 



SI = a+CLi- 1/2, (1.9) 



and 



Ka+} = Sij, (1.10) 
{ai,aj} = 0. 



The Hamiltonian then becomes quadratic in the a* 



H = J±J2 \( a t a i+i + a i a t+i) ~ 9HbH ^{ajai - ]-). (1.11) 



Operators a, can be expressed in a& by a Fourier transformation 

a,i = -)= y^exp(iA:-R ; )a fc . (1-12) 
ViV . 

NS(R l — Rj) = ^2exp(ik(Ri- Rj)) : (1.13) 

{a,, a,,} = ^fc-fc') (1.14) 

where Ri is the position vector of spin sites i and N is the total number of sites. The 
Hamiltonian is then diagonalized 

H = J2^ata k + ^ B HN, (1.15) 
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and it describes a free fermion gas. The dispersion relation e& is 

e& = J± cos k — qhbH. (1-16) 

Thus the energy band is a cosine band and g^sH is the chemical potential. For zero 
field, the Fermi level is at tj = 0. All states below the Fermi level are filled and all 
states above the Fermi level are empty as shown in Fig. 1.1 (a). A hole in the Fermi 
sea, called a spinon, corresponds to angular momentum |. The elementary excitations 
are partical-hole pairs, carrying angular momentum h. They are composite structures 
and their energy is the sum of the energy of a quasi-particle (spinon) of wave vector 
k\ and that of a quasi-hole of wave vector k 2 = q — k\. Their dispersion relation is 
defined as 

e 9 = efci+e* 2 , (1-17) 

which is doubly degenerate since the Fermi surface consists of two points, | and ^l. 
For each wave vector q, there is a continuum of magnons corresponding to spinons as 
shown in Fig. 1.1 (b). The lower limit of the continuum is constructed at k\ = tv/2, 
k 2 = q — 7r/2 and the upper limit at ki=k 2 = 3 Y L 

7T 7T 

Cmm(q) = Jl (cos(-) + cos(& - -)) = J ± sink, (1.18) 
tmax(q) = J±(cos(— h cos( - )) = 2 J ± sin -. (1.19) 

The magnetized state of a spin- 1/2 chain can be mapped onto a spinon band which 
is more than half filled as shown in Fig. 1.1 (c). The corresponding low energy 
excitation spectrum is shown in Fig. 1.1 (d). It is again a continuum but with more 
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k/Jt q/Jt 

Figure 1.1: Exact zero and finite field results for the XY spin- 1/2 chain obtained 
through the Jordan- Wigner transformation. (a)Band structure for spinons in zero 
field. (b)Boundaries of the resulting two spinon continuum, (c) Band structure for 
spinons in a magnetic field, (d) Boundaries of the resulting spectrum of longitudinal 
excitations. The transverse modes are omitted for the sake of clarity. 
The Heisenberg spin chain has qualitatively similar properties and is discussed in 
detail in Ref. [12]. 



complicated kinematical limits shown in solid lines. Interestingly, spin chains develop 
soft modes at the incommensurate wave vector 

qi = nir ± 2ir(S z ) (1.20) 

where (S z ) is the magnetization. 

The Heisenberg S=l/2 chain, though it is more complex than an anisotropic XY 
chain, has qualitatively similar properties [12]. Table 1.1 contains a list of quasi-one 
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dimensional Heisenberg spin 1/2 systems. To study the striking effects of a magnetic 
field on a S=l/2 chain experimentally, it is necessary to select a model system that 
is highly one dimensional with an ordering temperature that is small compared to 
exchange interactions. In addition, the reduced field, g/j, B H/J, should be large for a 
H 10 T magnetic field, which is currently available for neutron experiments. At the 
same time, it is important that the low energy part of the excitation spectrum can be 
resolved with sufficient energy resolution using a cold neutron triple axis spectrometer. 
These conflicting requirements make it difficult to find a suitable material. Field 
induced incommensurate soft modes have been observed in Copper Benzoate just 
recently [37]. An unexpected field induced gap has also appeared in Copper Benzoate. 
The gap in the excitation spectrum of Copper Benzoate was believed to result from 
a staggered field due to an alternating gyromagnetic tensor and the Dzyaloshinskii- 
Moriya interaction. Based on this analysis, a quantum Sine-Gordon theory has been 
developed[47, 48, 49, 50, 51]. The Sine-Gordon model is exactly solvable and many 
physically relevant quantities can be calculated. Unfortunately, it is very difficult 
to grow large single crystals of Copper Benzoate to adequately study details of its 
excitation spectrum experimentally. In Chapter 3, we present heat capacity and 
elastic neutron scattering measurement on a quasi one dimensional S=l/2 Heisenberg 
antiferromagnet CuC^ • 2DMSO (CDC). The results demonstrate that CDC belongs 
to the same classification as Copper Benzoate and is a clean model system that is 
particularly well suited for magnetic field studies. 
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Chemical 


Common 


J 


T N 


9HbH/J 


Ref. 


Formula 


Name 


(meV) 


(K) 






Sr 2 Cu0 3 




189(17) 


5 


0.006 


[14]-[18] 


SrCuC-2 




181(17) 


<1.5 


0.006 


[18, 19] 


Ca 2 Cu0 3 




«86 


8 


0.013 


[20, 21] 


BaCu 2 Ge 2 07 




46.5 


8.5 


0.024 


[22, 23] 


BaCu 2 Si 2 7 




24.1 


9.2 


0.048 


[22], [24]- [26] 


KCuF 3 




17.5 


39 


0.067 


[27]-[31] 


CuGe0 3 




10.4 


T SP = 14 




[32]- [34] 


CuCl 2 • 2N(C 6 H 6 ) 


CPC 


2.31 


1.14 


0.50 


[35, 36] 


Cu(C 6 H B COO) 2 -3(H 2 0) 


Cu-Benzoate 


1.57 


0.8 


0.73 


[37, 38] 


CuCl 2 -2[(CH 3 ) 2 SO] 


CDC 


1.43 


0.91 


0.81 


[39, 40] 


Cu(C 4 H 4 N 2 )(NO s ) 2 


CuPzN 


0.90(1) 


<0.05 


1.29 


[41] 


CuS0 4 • 5(H 2 0) 




0.25 


0.100 


4.63 


[42, 43] 


CuSe0 4 • 5(H 2 0) 




0.15 


0.125 


7.72 


[42, 43] 



Table 1.1: Key characteristics of S=l/2 quasi-one-dimensional Heisenberg antiferro- 
magnets. Table is taken from Ref. [44]. J refers to the intra-chain exchange constant. 
T N corresponds to the temperature below which long range magnetic order develops. 
T S p indicates the critical temperature for a spin-Peierls transition. g(i B H/ J is calcu- 
lated with g = 2 and H = 10 Tesla. Further tabulations of S=l/2 HAFM chains can 
be found in references [45] and [46]. 
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1.3 Integer spin chains 

The spin Hamiltonian for an integer spin chain can be written as 

H = J S* • S,+i + ]>>(S/) 2 + E[(Sn 2 - (Si y ) 2 )] (1.21) 

i i 

In Eq. 1.21, the first isotropic term is the nearest neighbor Heisenberg exchange 
interaction. The second and third term are anisotropy terms. D and E are the 
out-of-plane and in-plane anisotropy parameters, respectively. 

The spin Hamiltonian in Eq. 1.21 is not exactly solvable. By mapping the Hamil- 
tonian to the 0(3) non-linear sigma model in the limit of large spins, Haldane showed 
that integer spin chains have properties that are fundamentally different from half 
integer spin chains [13]. While the latter have a gapless two-spinon continuum of 
excited states, integer spin chains have a resonant triplet excitation separated from 
a singlet ground state by a gap. The ground state is short range ordered even at 
zero temperature, with a correlation function that decreases exponentially at large 
distance [13, 52]: 

(S -S,) « (-l)'|/|- 1/2 exp(-|/|/£). (1.22) 

A simple description of the ground state of the antiferromagnetic Heisenberg spin- 
1 chain was proposed by I. Affleck, T. Kennedy, E. H. Lieb and H. Tasaki[53]. The 
ground state is a valence bond solid. Each spin-1 may be viewed as two spin- 1/2 
spins that combine on site in a fully symmetric triplet wave function. Each of these 
two spin- 1/2 spins then form an antisymmetric singlet wave function with a spin- 1/2 
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<r "3> <H "g> <H ~S> O ""5> 

(a) 




(b) 

Figure 1.2: (a) A schematic representation of valence-bond ground state for S=l 
chains, (b) Excited states of the spin-1 chain correspond to breaking a valence bond 
and creating a triplet in the valence-bond ground state. 

from either its right or left nearest neighbor. (see Fig. 1.2 (a)) The unique singlet wave 
function of the entire chain resulting from this construction is a non-magnetic singlet, 
as it is built from many singlet pairs. Excited states of the spin-1 chain correspond 
to breaking a valence bond in the valence bond solid state, (see Fig. 1.2 (b)) Broken 
valence bonds can propagate coherently through the chain where they only scattering 
from other thermally excited triplet "defects" in the valence bond solid. Creation of a 
broken valence bond at rest requires an energy of A = 0.41050(2) J for the bi-linear 
spin-1 chain. 

In the case of anisotropic chains for integer spin with fl / 0, the first excited 
state, which is a degenerate triplet is split into a singlet and a doublet. The gap 
energy was accurately determined [55]: 

A z = Ao + 1.4LD 

A xy = A -0.57D. (1.23) 
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A wealth of experimental and numerical studies have provided spectacular con- 
firmation of the Haldane conjecture and they have greatly advanced understand- 
ing of the physics involved[52, 54, 55, 56, 57, 58, 59, 60, 61]. Table 1.2 lists some 
of experimental model systems for the S=l antiferromagnetic spin chains. One of 
the few remaining unresolved issues is the behavior of Haldane-gap antiferromag- 
nets in high magnetic fields. An external magnetic field splits the triplet of Haldane 
excitations[72, 82, 83, 84, 85, 86], driving one of the modes to zero energy at a crit- 
ical field H c = A/g(i B . Above H c , the pure 1D-HAF system is expected to recover 
a gapless ground state at T=0, without long range magnetic ordering but with an 
infinite correlation length, as for the magnetized state of spin-1/2 chain. One of the 
goals of the present research has been to explore the possibility of incommensurate 
spin correlations in this magnetized phase of the spin-1 chain. 

In Chapter 4, I present our neutron scattering study of the S=l chain, 
Ni(C 5 Di 4 N2)2N 3 (PF 6 ) (NDMAP). We first performed an inelastic neutron scattering 
experiment at zero field to characterize the spin Hamiltonian for NDMAP. In an ap- 
plied magnetic field, we found commensurate transverse magnetic ordering by means 
of neutron diffraction. Depending on the direction of the applied field, the high field 
phase has either three dimensional long range order or quasi-two-dimensional short 
range order. We then performed extensive neutron scattering experiments in the 
magnetized state of NDMAP. While the experiments do not provide direct evidence 
for incommensurate correlations, the ^-dependence of quasi-elastic scattering in the 
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Chemical 


Common 


J 


\J'/J\ 


D/J 


T N 


gUB H/J 


Ref. 


Formula 


Name 


meV 


10" 3 


10" 3 


K 


10" 2 




AgVP 2 S 6 




58(4) 


0.01 


5.8 


< 2 


2 


[62] 


Y 2 BaNi0 5 




21 


< 0.5 


-39 


< 0.05 


5.5 


[63]- [65] 


Ni(C3Hi0N 2 )2N 3 (ClO4) 


NINAZ 


10.7 


< 0.7 


170 


< 0.06 


11 


[66]- [68] 


Ni(C 3 Hi0N2)2NO 2 (ClO4) 


NINO 


4.5 




210 


< 1.2 


26 


[69, 70] 


Ni(C 2 H 8 N2)2N0 2 (C104) 


NENP 


4.1(3) 


0.8 


180 


< 0.0003 


28 


[71]-[73] 


Ni(C 6 Di4N2)2N 3 (PF 6 ) 


NDMAP 


2.85 


0.6 


250 


< 0.25 


41 


[74]- [78] 


CsNiCls 




2.275 


17 


-1.9 


4.9 


45 


[79]-[81] 



Table 1.2: Key characteristics of quasi-one-dimensional an S = 1 Haldane antifer- 
romagnets. J refers to the intra-chain exchange constant. J' refers the inter-chain 
exchange constant. D corresponds to the single ion anisotropy parameter. T N refers 
to the temperature below which long range magnetic order develops. g/i B H/J is 
calculated with g = 2 and H = 10 Tesla. 



high field phase is consistent with such a possibility. 

In Chapter 5, I will discuss impurity effects in the S=l quasi-one dimensional 
Heisenberg antiferromagnetic chain Y 2 BaNi0 5 . One scenario is to replace the mag- 
netic Ni 2+ ions by non-magnetic ions such as Zn 2+ or Mg 2+ . Since those impurities do 
not have magnetic moments, one would not expect interactions between Ni 2+ ions to 
propagate through these impurity sites. This type of chain-cutting effect was observed 
by DiTusa et aim their studies of powder samples of Zn doped Y 2 BaNi0 5 [64], and by 
Xu et al in their studies of Mg doped Y 2 BaNi0 5 [87]. Another scenario is to replace 
the off-chain Y 3+ by Ca 2+ , which gives rise to far more interesting cooperative effects. 
Ca-doping introduces holes onto the chains although the material does not become 
metallic. Inelastic neutron scattering experiments reveal the existence of excitations 
below the gap with an incommensurate structure factor[88] that reflects the internal 
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structure of the spin polarons that develop around the hole. We present our inelastic 
neutron scattering study of excitations in the bulk of the S=l chain and excitation 
associated with spin polarons in quantum magnets Y 2 _ ;r Ca ;r BaNi05 (£=0.04, 0.1). 
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Chapter 2 

Experimental techniques 

2.1 Crystal Growth 

For most neutron scattering experiments, single crystal samples are needed in 
order to fully investigate the dynamic properties of the systems. 

Our collaborators from Japan, Dr. Koichi Katsumata and Dr. Zentaro Honda, 
provided deuterated single crystal samples of NDMAP investigated in Chapter 4. 
Dr. Ito provides single crystal samples of Y 2 BaNi0 5 and Ca doped Y 2 _ x Ca x BaNi0 5 
(x=0.04, x=0.1) investigated in Chapter 5. I grew deuterated CDC samples used in 
the experiments described in Chapter 3. 

Deuterated single crystals of CuCl 2 • 2((CD 3 ) 2 SO) (CDC) were obtained by slow 
cooling from the temperature T = 50 °C to T = 20 °C of saturated methanol solutions 
of anhydrous CuCl 2 and deuterated dimethyl sulfoxide in 1:2 molar ratio. Here we 
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describe one of the procedures that we followed to obtain the crystals. The chemical 
equation that we used is 

CuCl 2 + 2((CH 3 ) 2 SO) -> CuCl 2 • 2((CH 3 ) 2 SO). (2.1) 

We dissolved 2.68 g of CuCl 2 into 30 ml of methanol at 50°C, then added 2.85 ml 
of dimethyl sulfoxide. It is then well mixed to make a bright, blue solution. The 
solution was held at 50° C for an hour. Over the next 10 hours, the temperature was 
reduced to 40° C. Over the next 100 hours, the temperature was lowered to room 
temperature. Around 4 g of Emerald-green CDC crystals grow on the bottom of the 
beaker as large tabular plates with well developed (001) faces. 



2.2 Heat Capacity 

The specific heat of a system is defined as 

C - (™) 
~ dT 

= — J en(e)D(e)de, (2.2) 

where E is the internal energy, n(e) is a thermal population factor, and D(e)) is the 
density of states. By measuring low temperature specific heat, one obtains informa- 
tion about the low energy lattice and electronic excitation spectrum. Low temper- 
ature specific heat measurements are also useful for revealing phase transitions in 
magnetic and superconducting materials. For this dissertation, the low temperature 
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specific heat of a spin 1/2 quasi-one-dimensional antiferromagnet was measured to 
study magnetic phase transitions and the field induced spin excitation gap of the 
system. 

Low temperature specific heat was measured by the time relaxation method, which 
was developed in 1972 by Bachmann et al. [89]. This method is particularly well suited 
to measure small samples (~ 10 mg) at low temperature (0.1 K < T < 16 K). The 
basic idea of the time relaxation method can be understood by an electrical analogue 
as illustrated in Fig. 2.1. The sample is mounted on a calorimeter and is connected 
to the heat reservoir by a weak thermal link with thermal conductance k as shown 
in Fig. 2.1 (a). The specific heat of the sample and calorimeter, its temperature, 
and 1/k can be mapped onto the capacitance, c, the voltage across the capacitor, 
and the resistance R of the parallel resistor-capacitor circuit shown in Fig. 2. 1(b). In 
steady state, power, P, flows through the calorimeter and this corresponds to current 
/ flowing through the analogous electrical circuit. When we turn off the power, the 
temperature T exponentially relaxes to T : 

T = T + ATe- (< - <o)/ri (2.3) 

with a characteristic time constant 

n = f (2.4) 

where AT = T — T . The specific heat, C, is therefore determined by C = kt\. 

To measure specific heat, a step function in heater power is used. First the power 
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sample 




V = iR 



calorimeter Q 'J 




1/k 



c 



R 



V = 



thermal 



-L 



(a) 



(b) 



Figure 2.1: Thermal relaxation method: (a) schematic heat flow circuit (b)electrical 



is adjusted such that the sample plus calorimeter is heated a small AT above the 
thermal bath temperature T (AT/T ~ 1%). Then the power is turned off, and 
the temperature of sample plus calorimeter decays exponentially. The value of the 
thermal conductivity k is given by Pj AT. t\ is determined by fitting the measured 
exponential curve to Eq. 2.3. The heat capacity of sample plus calorimeter then 
is determined by C = kt\. By subtracting the heat capacity of the calorimeter, we 
obtain the heat capacity of the sample. To improve signal to noise ratio, several power 
pulses are applied, and the average of signals from each pulse is used to determine k 
and T\. 

In the case that the thermal conductance between sample and the calorimeter or 



analog 
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the internal thermal conductance of sample is poor, a second time constant r 2 should 
be taken into consideration [90]. 

2.3 Neutron Scattering 

2.3.1 Introduction 

The neutron is an uncharged particle with mass 1,839 times that of the electron. 
The life time of a free neutron is around 998 seconds, although it is a stable particle 
when bound in an atomic nucleus. The neutron is a spin-1/2 particle and therefore 
has a magnetic dipole moment 1.913/iijv (//jv=5.0508 x 10~ 27 J/T). These unique 
properties of the neutron make it an ideal probe of condensed matter systems. 

The zero net charge of the neutron means that neutrons interact directly with 
nuclei rather than with the electronic cloud. This has three important consequences: 
1) it is possible to study light elements, in the presence of heavier ones: 2) neutrons 
can easily distinguish atoms with comparable atomic number; 3) neutrons distinguish 
isotopes. The interaction of a neutron with the nucleus of an atom is rather weak, 
making it a highly penetrating probe. This also allows neutron to be a non-destructive 
probe of the interior of materials. 

Cold and thermal neutrons have energies ranging from a few meV to a few hundred 
meV, which is similar to the energies of atomic and electronic processes in condensed 
matter. This makes neutron ideal for studying a wide variety of dynamical phenomena 
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in solids and liquids. At the same time, the de Broglie wavelength of cold and thermal 
neutrons is comparable to the atomic spacing in solids so that interference phenomena 
that can reveal spatial structure and correlations are also possible. 

The magnetic moment of the neutron makes it very useful for probing magnetic 
structure and dynamics of a magnetic system. The cross-sections for magnetic and 
nuclear scattering are of the same magnitude, making it possible to simultaneously 
measure the magnetic and chemical structure. The 1994 Nobel Prize in Physics 
was awarded to Bertram Brockhouse and Clifford Shull for their pioneering work in 
developing neutron scattering techniques for studying condensed matter, in particular, 
for the development of the neutron diffraction and neutron spectroscopy techniques. 

2.3.2 the Scattering of Neutrons 

All neutron scattering experiments are essentially measurements of the partial 
scattering cross section, d 2 a / dfldE f , as a function of neutron momentum transfer. 
Q, and energy transfer, huj. A neutron incident on a sample with wave vector kj and 
energy Ei is scattered into a final wave vector k/ and final energy Ef. Because of 
conservation of momentum and energy, the momentum 

Q = k i -k / (2.5) 

and the energy 
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Q = k f -k f 

20 > ha = E, - E f 

Figure 2.2: Graphical presentation of the momentum conservation law. \Q\ = k\ + 
cos 29 where 29 is the scattering angle. Such a diagram is also known as a 
scattering triangle. 



} 2m 2m y ' 

are transferee! to the sample. This can be represented by the scattering diagram 
shown in Fig. 2.2. 

Magnetic neutron scattering occurs because of the electromagnetic dipole-dipole 
interaction between neutrons and unpaired electrons in the atom. The partial differ- 
ential scattering cross-section describing the scattering from a Bravais crystal with 
localized spins by unpolarized neutrons can be written as [91, 92] 

-J^E = (7ro)2 ^f F (Q)\ 2 e- 2W(Q) " QMS^{CtM (2.7) 

f 1 a/3 

where 7 = 1.913 is gyromagnetic ratio of the neutron, r is the classical electron radius 
(e 2 /m e c 2 = 2.82 fm), e~ 2W ^ is the Debye- Waller factor, -F(Q) is the magnetic form 
factor [93]. S a P(Q,uj) is magnetic scattering function, which corresponds to the 
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Fourier transformed two point dynamics spin correlation function: 



S^(Q,co) = Y*[ dte-^j-^e^^ < S?(0)Sl(t) > . (2.8) 
77 ^ i',i 

The magnetic scattering function is also directly related to the imaginary part of 
the generalized spin susceptibility via the fluctuation-dissipation theorem 



S(Q, U )= , ' (2.9) 

where f3 = ljksT. The fluctuation-dissipation theorem allows direct comparison of 
S a/3 (Q,uj) attained by neutron-scattering to \ obtained by bulk measurements or 
theoretical calculations. 

Replacing elements in Eq. 2.7 with their thermal averages as t — >■ oc and inte- 
grating with respect to Ef gives the elastic cross section. As t — > oo, < ^"(O)^^) > 
becomes independent of time and is equal to < >< Sy >. Substituting this in 
Eq. 2.7 and integrating with respect to Ef, the elastic cross-section is given by 



^ = (7ro) 2 ||F(Q)| 2 e- 2 ^ £(<W " QaQ/0£ e<Q ' (1 '" 1) < Sf >< S{ > . (2.10) 

a/3 l',l 

Elastic scattering from a magnetical crystal in which the spins are ordered gives rise 
to magnetic Bragg peaks and is used to determine the orientation and magnitude of 
the atomic magnetic moments, i.e., the magnetic structure of a crystal. In this case, 
Eq. 2.10 becomes 
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% = (jro) 2 N m { ^- J](|F(Q)| 2 - |Q ■ F(Q)| 2 )5(Q - r m ). (2.11) 
The magnetic vector structure factor is defined as 

F(Q) = ]T f F d (Q)e- 2W ^ <S d > e^ d (2.12) 

d 

In the above two equations, N m is the number of magnetic unit cell, v m is the 
volume of magnetic unit cell, r m is a vector in the magnetic reciprocal lattice, d the 
position vector of each ion within the magnetic cell. i*d(Q) is the magnetic form 
factor for ion at position d. 

2.3.3 Instrumentation 

Measurement of the neutron scattering cross section requires controlled access to 
wave vector transfer Q and energy transfer hco, which can be achieved in various ways. 
The first and simplest type of instrument is the triple-axis spectrometer. The three 
axes correspond to rotation of the monochromator, the sample, and the analyzer, as 
shown in Fig. 2.3. The monochromator crystal selects neutrons with specific incident 
wavelength from a "white" neutron beam, to arrive at the sample through Bragg 
diffraction 



= 2d M sin 6 M 



(2.13) 
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where n is a positive integer, d m is the crystallographic d- spacing for the monochro- 
mator crystal. During the interaction with the sample, neutrons may lose or gain 
energy, so they emerge with an energy different from the incident energy. The ana- 
lyzer selects another specific energy for scattered neutrons to be detected, again via 
Bragg's law 

nX f = 2d A sm9 A . (2.14) 

The material most commonly used as monochromator and analyzer in a triple-axis 
spectrometer is pyrolytic graphite (PG). This is due to its high reflectivity and broad 
mosaic distribution for neutrons over a wide range of energy. Beryllium, Copper, 
Germanium, and Zinc are used as monochromator crystals for higher energy neutrons. 

We usually use the neutrons with the fundamental (n=l) wavelength A/. However 
neutrons with higher-order wave length (A = Xf/n,n > 1) also satisfy the Bragg 
condition and will also be present in the Bragg diffracted beam. These energetic 
neutrons can pose a serious problem as background contamination. Filters are placed 
in the path of neutron beam to remove these unwanted higher order neutrons. PG 
is widely used as a filter for thermal neutrons. For cold neutrons, polycrystalline 
beryllium is commonly used. Both filters function by Bragg diffracting neutrons 
away from the direction of transmission. 

The triple-axis-spectrometers used in experiments presented in this dissertation 
are the BT2, and BT4 thermal neutron spectrometer and the cold spectrometer 
SPINS shown in Fig. 2.4. All these instruments are at the NIST center for neutron 
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Figure 2.3: Schematics of a three-axis spectrometer. Figure is taken from R. M. 
White, Quantum Theory of Magnetism, Springer (1970). 
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Figure 2.4: The SPINS cold neutron spectrometer at the NIST Center for Neu- 
tron Research in Gaithersburg, MD, USA. The blue cylindrical object marked Ox- 
ford, is a high-field dilution refrigerator. The figure is taken from NIST webpage 
http: / / ncnr.nist.gov / instruments / spins/. 



research at the National Institute of Standard Technology, Gaithersburg, Maryland. 

Triple-axis-spectrometers are very efficient for probing a few points in the (Q, u) 
space. If data are required for a large (Q,w) range, a time-of-flight spectrometer is 
used. The advantage of this technique over the crystal spectrometer is that large 
regions of energy transfer and momentum transfer can be measured simultaneously 
with the use of large arrays of position-sensitive detectors, resulting in a gain in 
efficiency. In the direct geometry time-of-flight method, a pulsed monochromated 
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neutron beam with energy Ei incident on the sample. The monochromatic neutron 
pulses are then scattered from the sample and detected by a large array of detectors. 
Each of these detectors is capable of collecting a time dependent spectrum so that the 
flight time, t, can be determined for each detected neutron over flight path of length 
L. The final energy of the detected neutron is determined by Ef = \mv^ = \m{ j) 2 . 
Thus the energy transfer huj = Ei — Ef is determined. The modulus of the wave 
vector transfer Q is determined by |Q| = ^jk\ + k 2 — 2kikf cos 29, where 26 the 
scattering angle at which the detector is positioned. The vector components of Q can 
be determined by the azimuthal angle of the detector. 

The time-of-flight instruments used for this dissertation was the MAPS Fermi 
chopper spectrometer at the ISIS facility, Rutherford Appleton Laboratory, United 
Kingdom, which is shown in Fig. 2.5. 
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Figure 2.5: Schematic diagram of the MAPS Chopper Spectrometer at the ISIS fa- 
cility Rutherford Appleton Laboratory, United Kingdom. One of the innovations in 
MAPS is a 16 m 2 position sensitive detector array which provides 160,000 detector 
pixels with almost continuous coverage in the forward angle direction. This figure is 
taken from the ISIS webpage at http://www.isis.rl.ac.uk. 
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Chapter 3 

S=l/2 Heisenberg 
antiferromagnetic chain: CDC 

3.1 Introduction 

The dynamical behavior of the Heisenberg antiferromagnetic spin- 1/2 chain has 
been extensively studied. The excitation spectrum is a bounded continuum resulting 
from excitation of fermionic spinons. Field dependent soft modes have been pre- 
dicted [12, 94, 95] and observed in a recent neutron scattering experiment on Copper 
Benzoate in high magnetic field [37]. However, an unexpected gap in the excitation 
spectrum was induced by the applied field. The gap was theorized to result from a 
staggered field due to an alternating gyromagnetic tensor and the Dzyaloshinskii- 
Moriya interaction. Based on this analysis, a quantum sine-Gordon theory was 
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developed [47, 48, 49, 50, 51]. The quantum sine-Gordon model has successfully ex- 
plained the specific heat, susceptibility and field induced gap of copper benzoate, 
[PM • Cu(N0 3 ) 2 • (H 2 0) 2 ] n [96] and Yb 4 As 3 [97]. The sine-Gordon model is exactly 
solvable and many physically relevant quantities can be calculated. In particular, 
the spectrum is known to consist of a soliton-antisoliton massive doublet and their 
bound states, which are called " breathers". It is very difficult to grow large single 
crystals of copper benzoate to adequately study details of its excitation spectrum 
experimentally. In this work, we present heat capacity and elastic neutron scattering 
measurement on CuCl2 • 2DMSO (CDC). The crystal structure of CDC is similar to 
that of Copper Benzoate in that it consists of one dimensional chains of S=l/2 Cu 2+ 
ions that have an alternating g-tensor along the chain axis. However, large single 
crystals of CDC can readily be grown from solution. Here I show through heat ca- 
pacity and elastic neutron scattering study that CDC belongs to the same class as 
Copper Benzoate and is a clean model system that is particularly well suited for an 
experimental study of the quantum sine-Gordon model. 

3.2 CDC 

CuCl 2 • 2((CD 3 ) 2 SO) (CDC) was previously identified as a spin-1/2, linear antifer- 
romagnetic chain with exchange constant J = 1.43 meV[40]. The structure of CDC is 
orthorhombic, with space group Pnma, and lattice constants a = 8.054 A, b = 11.546 
A, and c = 11.367 A[98]. The Cu 2+ ions form chains along the a axis as shown in 
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Fig. 3.1 . The local symmetry axis of the copper sites lies in the a-c plane, but is 
tilted by ±22° from the a axis for alternate sites. This alternation of local crystalline 
environment leads to an alternating g tensor in CDC. We assume that g tensor takes 

the form[40] similar to what found in copper benzoate[99]: 

/ \ 

2.28 ±g s 



1.97 



V 



g u ±g s , 



(3.1) 



±9* 2.12 

where g u and g s are the uniform and staggered parts of the g tensor. The spin 
Hamiltonian describing a Heisenberg in S=l/2 chain in magnetic field applied along 
z direction can be written as 



H = J^Si-S i+1 - wFLS*. 



(3.2) 



However, in CDC, the alternating g tensor results in the presence of an effective 
staggered field, upon the application of a uniform field. The second term in Equation 



3.2 must be replaced by 



mag 



i,a,b 



(3.3) 



At the same time, a staggered Dzyaloshinskii-Moriya (DM) interaction also con- 
tributes to the effective staggered field[47, 48, 100, 101]. 



H DM = J2(-iyn-(Si-i xso 



(3.4) 



The direction of the DM vector, D, can be constrained by considering the symmetries 
of the crystal structure of CDC. First, the structure is invariant under a translation 
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along the a axis by two Cu sites. This means that there can be at most two dis- 
tinct DM vectors. Second, the a-c plane is a mirror plane, which implies that the 
DM vectors point along b[101]. Finally, the crystal structure is invariant under the 
combined operation of translation by one Cu site along the chain, and reflection in 
the a-b plane. Thus, the DM vectors alternate along the chain as D = ±Db, and the 
contribution to the spin Hamiltonian from the DM interaction is 

n DM = • (Si_! x S^. (3.5) 

i 

Thus both the alternating g tensor and the alternating DM interaction give a stag- 
gered field h when a uniform external field is applied. In general, h is nearly orthog- 
onal to H and can be approximately given by 

h^Dxs'Hlj'H, (3.6) 

Taking into account that h « H and h is almost perpendicular to H, the spin 
Hamiltonian of CDC in a magnetic field can be written as 

# = ^ JS ; • Si+i - WHS? - iMsi-lYhSf, (3.7) 

i 

which corresponds to a Heisenberg model with orthogonal uniform and staggered 
fields. 

3.3 Experimental procedure 

Single crystals of CDC were obtained by slow cooling from temperature T = 50 
°C to T = 20 °C of saturated methanol solutions of anhydrous CuCl 2 and dimethyl 
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sulfoxide in 1:2 molar ratio. Emerald-green crystals grow as large tabular plates with 
the (001) faces well developed. (See Section 2.1 for details of crystal growth.) Spe- 
cific heat measurements in both the ordered and disordered phases were performed on 
small, protonated crystals with a typical mass of 15 mg in a dilution refrigerator using 
relaxation calorimetry in magnetic fields H up to 9 T. To determine the magnetic 
structure of the ordered phase, neutron diffraction experiments were performed on 
deuterated single crystals of mass ~ 1 g using the SPINS cold-neutron triple-axis spec- 
trometer at the Center for Neutron Research at the National Institute of Standards 
and Technology. The spectrometer was configured with horizontal beam collimations 
of 50' /ki - 80' - 80' - 300'. A pyrolytic graphite (PG(002)) monochromator was used 
to select incident neutron energies of either E{ = 3 or 5 meV, and a cooled Be filter 
was used before the sample to eliminate contamination from higher-order reflections. 
Measurements were performed in both the (hkO) and (hOl) scattering planes at tem- 
peratures T > 0.25 K in magnetic field up to H = 7.25 T applied perpendicular to the 
scattering plane. The low-temperature magnetic structure was determined from rock- 
ing scans of magnetic Bragg peaks. An absolute scale for the integrated intensities 
of these reflections was obtained from measurements of nuclear Bragg peaks, and the 
the magnetic structure factors squared were obtained after correcting for resolution 
effects, the sample mosaic spread and the magnetic form factor of the Cu 2+ ions. 
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Figure 3.1: Crystal structure of CDC. ac plane and a three dimensional cell are shown. 
a axis is the chain axis. The axes of the quasi-planar CuCl 2 2 groups alternate 22° 
on either side of the a axis. 
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3.4 Experimental results and data analysis 



3.4.1 Specific heat 

Fig. 3.2 shows the temperature dependence of the total specific heat of CDC for 
a number of magnetic fields applied along the crystalline a, b, and c axes. In zero 
field, the peak at 0.91 K corresponds to three-dimensional long range ordering. Fig. 
3.3 shows the phase boundaries derived from these data and neutron diffraction data 
that will be discussed later. When the field is applied in the mirror plane, either 
along the a or the c axis, the Neel temperature decreases as the field increases. When 
the field is increased above a critical value, H c m 3.9 T for fields applied along c axis, 
and H c m 5.7 T for fields applied along the chain axis, the specific heat becomes 
exponentially activated with spectral weight shifting to higher temperature region as 
the field is further increased. This clearly indicates a field induced gap A(iJ). When 
the field is applied perpendicular to the mirror plane, i.e. along the b axis, however, 
the Neel temperature increases as the field is increased. We fit the H > H c data in Fig. 
3.2 (a) and (c) to explore field dependence of the energy gap. The low temperature 
specific heat is first fitted to a simple model of h species of one dimensional bosons 
with a gap as shown in Fig. 3.4[37, 102]. The extracted gap values for H\\a and 
H\\c are shown in Fig. 3.5, the solid lines are the power law fits. The specific heat 
data is further fitted to a quantum sine-Gordon model shown in Figure 3.6. Figure 
3.7 shows energy gap versus field for H\\a and H\\c respectively. Solid lines are the 



36 

theoretical result for sine-Gordon model that don't fit our data well for reasons that 
will be explained later. 

3.4.2 Neutron diffraction 

The (h,0,l) and (h,k,0) reciprocal lattice zones were searched for magnetic Bragg 
peaks. Fig. 3.8 shows a map of the reciprocal space that we covered in our experiment. 
As a result, magnetic Bragg reflections of (h,k,0) with h odd integers and k integers 
were observed at temperature below the Neel temperature. All wavevector transfers 
in terms of reciprocal-lattice unit of the crystallographic structure. Only the (1, 0, 0) 
and (3, 0, 0) magnetic reflections were observed in the (h,0,l) zone. Fig. 3.9 shows (h, 
0, 0) scans through (1, 0, 0) magnetic peak under different conditions of temperature 
and applied magnetic field. At zero field and T = 0.3K, below T N = 0.91K, a 
well defined resolution limited peak is observed centered at h = 1.0. When the 
temperature is increased above the Neel temperature, T = 2K, the peak vanishes, 
thus clearly indicating that it represents the order parameter of the transition. When 
a field is applied along the b axis at T = 0.25K, the magnetic intensity at (1, 0, 
0) is enhanced compared to the zero field data. The extra intensity at h ~ 1.1 is 
temperature independent, therefore it can not be of magnetic origin but rather a 
spurious effect. No incommensurate order was observed throughout the experiment. 

To extract the critical exponent (5 of the magnetic ordering transition, we mea- 
sured the temperature dependence and field dependence of magnetic Bragg peak 
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Figure 3.2: Specific heat of CDC as a function of temperature for several magnetic 
field values. The magnetic field is applied parallel to the (a) a, (b) b, and (c) c axes, 
respectively. 
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Figure 3.3: Field-temperature phase diagram of CDC determined from heat capacity 
measurements (open symbols) and from the neutron scattering measurements (filled 
symbols). The diamonds represent a spin-flop phase transition for H\\c as observed 
in neutron scattering experiment. The dashed lines are a guide to the eye. 
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Figure 3.4: Semilogarithmic plot of C(H,T) vs 1/T in the gapped phase with fits 
to non-interacting a one-dimensional boson model as described in the text. The 
magnetic field is applied parallel to the (a) a and (b) c axes, respectively. 
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Figure 3.5: Field dependence of the energy gap derived from fits to specific heat data 
such as shown in Figure 3.4. The solid lines are from fits to power laws described in 
the text. 
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Figure 3.6: Semilogarithmic plot of C (H,T) vs 1/T in the gapped phase with fits to 
sine-Gordon model as described in the text. The magnetic field is applied parallel to 
the (a) a and (b) c axes, respectively. 



42 




Figure 3.7: Field dependence of the energy gap derived from fits to specific heat data 
such as shown in Figure 3.6. The solid lines are from fits to power laws described in 
the text. 
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Figure 3.8: (h, 0, 1) and (h, k, 0) reciprocal lattice zones were searched for magnetic 
Bragg peaks. The open and solid symbols represent the locations of the nuclear and 
magnetic Bragg peaks we observed, respectively. The dashed lines represent the limit 
of experimentally accessible portion of the reciprocal space. 
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Model Ginzburg-Landau 


2-d Ising 


3-d Ising 2-d X-Y 


3-d X-Y 


Heisenberg 


P 0.5 


0.125 


0.326 


0.345 


0.367 



Table 3.1: Approximate values of critical exponents ft for various models. See Mal- 
colm F. Collins, Magnetic Critical Scattering, Oxford University Press, (1989) for a 
complete list. 



intensity. We fit them to 

I(T) oc \T-T N \ 2 P (3.8) 
1(H) a \H-H C \ 2 ? 

respectively. Table 3.1 lists theoretical values of £ for various models. Figure 3.10 
shows the temperature dependence of the (100) magnetic reflection for different field 
magnitudes and orientations. The intensity decreases as the temperature approaches 
the Neel temperature T N . The solid line shows the result of a least-square fit to 
determine the critical exponent of £ and critical temperature T N . At zero field, 
we found (3=0. 301(1) and T N =0.91 K. The Neel temperature is consistent with the 
specific heat measurement. With a 7 T field applied along the b axis, T c = 1.07 K 
and £ = 0.351(1) were determined. When a 3 T field is applied along c axis, we 
get T c = 0.768K and £=0.363(2). Figure 3.11 shows field dependence of the (1, 0, 
0) magnetic reflection for different temperatures and field orientations. When the 
field is applied along the c axis, the intensity decreases as the field increases. A 
Power law fit gives the result of £=0.400(1) and # c =3.9 T. The step in intensity 
near 0.3 T is identified as a spin-flop transition. The phase digram for the spin-flop 
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transition was mapped out by field scans as a function of temperatures and it is 
shown by filled diamond symbols in Fig. 3.3. When a field is applied along b axis, 
the magnetic intensity of the (1, 0, 0) reflection is enhanced as shown in T = 0.25 K 
data. At T = 1 K, the field drives CDC to the ordered phase at H c ~ 4.05 T with 
f3 = 0.481(3), which is indistinguishable from mean field value j3 = 0.5, indicating 
the true critical region is too narrow to access here. Figure 3.12 shows logarithmic 
plot of the magnetic intensity (1, 0, 0) against reduced temperature or reduced field, 
(1— T/T c ) or \1—H/H C \, around the critical region. Almost all of the critical exponents 
(3 falls in the range of 0.30 ~ 0.36, and they are independent of the orientation of 
magnetic field, the magnitude of the applied field and the temperature as shown in 
Figure 3.12. The average value, /3=0.354, agrees well with /3=0.367 as determined for 
the 3D Heisenberg model. 

3.5 Discussion 

3.5.1 Specific heat 

The qualitative behavior of the specific heat is closely related to the crystal struc- 
ture of CDC. When there is an applied field in the ac plane, the alternating g tensor 
and the DM interaction lead to an effective staggered field. This staggered field re- 
sults in an energy gap. When the field is applied along the b axis, the two copper 
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0.8 0.9 1.0 1.1 1.2 

h (rlu) 

Figure 3.9: (h,0,0) scans through (1, 0, 0) magnetic Bragg peak in CDC at T = 0.25K, 
H = 7.25 T||6 axis, T = 0.25K and T = 2K at zero field. The bar in the figure 
represents the calculated peak width of (100). 



47 




Figure 3.10: Measured temperature dependence of the (1. 0. 0) magnetic peak in- 
tensity for zero field, H = 7.25 T applied parallel to the b axis and H = 3T applied 
parallel to the c axis. 
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Figure 3.11: Measured field dependence of the (1, 0, 0) magnetic peak intensity at 
T = 0.25 K with field applied along c axis and at T = 0.25 K and at T = 1 K 
with field applied along the b axis. The step in intensity near H = 0.3 T for H\\c 
corresponds to a spin-flop transition. 
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Figure 3.12: Logarithmic plot of the background corrected (1, 0, 0) magnetic intensity 
plotted as a function of reduced temperature (open symbols) or reduced field (solid 
symbols). 
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sites experience the same field and there is no staggered field in the crystal. We can 
also see this through Eq. 3.6. When H||6, D x g u H = g s H=0, thus staggered field 
h=0. 

To explore the details of field induced gap, the field-dependent specific heat was 
first analyzed in a naive model which was used in the analysis of Copper Benzoate 
[37]. We fit the high-field specific data in Fig. 3.4 to 

C ma9 = ^L(^-f>-exp(-A/k B T), (3.9) 
V27T k B T v 

which is the low-T specific heat [102] for h species of noninteracting one dimensional 
bosons with a gap and dispersion relation 

huj(q) = A + {vqf/2. (3.10) 

We fit the total specific heat of a system in all the fields to a global fit 

C(H) = C mag (H) + aT 3 + b^, (3.11) 

where aT 3 and bj^ are included to take into account the small lattice contribution 
and nuclear contribution[103], respectively. The fit results are shown in solid lines in 
Fig. 3.4. The fit gives a =0.0284 J/mol-K 4 , b=2.293 xlO~ 5 Jmol • K/T 2 . The values 
of v/n= 0.22(2). 0.18(3) meV/mode for H\\a and H\\c respectively, are close to the 
crude estimate v/n m O.bnJ/Q = 0.37 meV/mode, where n = 6 is the number of 
soft modes in the Brillouin Zone. These values are also close to the values reported 
in Copper Benzoate. The field dependence of the gaps derived from this analysis is 
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H||a 




H||c 


H (T) A (meV) 


H(T) 


A (meV) 


6 0.214(1) 


4 


0.250(5) 


7 0.268(5) 


5 


0.302(4) 


9 0.334(5) 


7 


0.409(6) 




8 


0.478(7) 




9 


0.505(4) 



Table 3.2: Field dependence of energy gap derived from a global fit of specific heat 
data to Eq. 3.9, which describes the low temperature specific heat for noninteracting 
one-dimensional bosons. 

summarized in Table. 3.2. and plotted in Fig. 3.5. The field dependence of the gap 
can be described by the power law A(H) = AH a with a = |, as found in Copper 
Benzoate. We find A a = 0.08(1)# 2 / 3 for H || a and A c = 0.11(1)# 2 / 3 for H || c if A 
is measured in meV and H in Tesla. 

Encouraged by the above results, we analyze our specific heat data in the frame- 
work of the quantum sine-Gordon model. The Hamiltonian described by Eq. 3.7 can 
be transformed by Abelian bosonization into a sine-Gordon model with Lagrangian 
density [47] 

L = ^(d u (i ) ) 2 + \(h)cos(pO) (3.12) 

where 6 is the dual field and the coupling f3 depends on the value of applied magnetic 
field. The spectrum of the sine-Gordon theory is well known: it consists of a soliton- 
antisoliton doublet of mass M. Their bound states are called breathers. The specific 
heat of the sine-Gordon model can be obtained by solving a set of integral equations 
based on Bethe Ansatz. It successfully explained the specific heat of Copper Benzoate 
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H||a 




H||c 


H (T) A (meV) 


H(T) 


A (meV) 


6 0.153(3) 


4 


0.206(4) 


7 0.186(4) 


5 


0.276(4) 


9 0.336(4) 


7 


0.367(4) 




8 


0.520(6) 




9 


0.640(5) 



Table 3.3: Field dependence of energy gap derived from a global fit of specific heat 
data to Eq. 3.13, which describes the low temperature specific heat of the quantum 
sine-Gordon model. 

under an applied field[50]. At low temperature the specific heat is found to be of the 
form 

°- ~E^[ 1 + i§7 + i^O* «p<-*./*.r>. (3.i3) 

o 1 v s 

M a is the mass of breathers and can be expressed in terms of the soliton mass mass 
gap M by 

M a = 2Msin(n</2) ; (3.14) 
n = l,....,[l/£], 

where £ = g^rrp-, and [l/£] denotes the largest integer smaller than l/£. The soliton 
mass, M, is [48] 

M/J~ 1.504 A 2 / 3 | In 4| 1/6 - (3.15) 



The field-dependent specific heat data is fitted to Eqs. 3.13 and 3.11 in a global fit, 
shown in Fig. 3.6 as solid lines. It gives a=0.0355 J/mol-K 4 , b=8.9 x 10~ 6 mol ■ K/T 2 . 
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Apparently, Eq. 3.13 describes the observed temperature dependence of the specific 
fit very well. A decrease of spin wave velocity by the application of the external field 
was observed as shown in Fig. 3.13. We get v s (6 T) : v s (7 T) : v s (9 T) = l : 0.97 
: 0.42 for H\\a and v s (4 T) : v s (5 T) : v s (7 T) : v s (8T):v s (9 T) = 1 : 0.70 
: 0.51 : 0.25 : 0.17. The decrease is faster than the theoretical predication, which 
gives 1 : 0.98 : 0.95 and 1 : 0.99 : 0.98 : 0.96 : 0.95 respectively [48]. Table 3.3 
and Fig. 3.7 summarize the field dependence of the gaps. At lower fields, the gap 
values are similar to those obtained by fitting the specific heat data to noninteracting 
one-dimensional bosons. However, at higher fields (H>8 T), the former values are 
consistently higher than the latter. The theoretical field dependence of soliton mass 
of sine-Gordon model is described by Eq. 3.15. Taking into account that the effective 
staggered field h = cH and the logarithmic term in Equation 3.15 varies only slowly 
with H, the field dependence of the soliton mass can be simplified as 



The solid lines in Fig. 3.7 are the power law fits to the gap values. The fits give 



M(H)/J = cH 1 * = 1.504(c/J)i#I. 



(3.16) 



M = 0.05(1)# 2/3 



(3.17) 




M = 0.10(1)# 2/3 



(3.18) 
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Figure 3.13: The spin wave velocity as extracted from fitting the specific heat data 
to Eq. 3.13. The data are plotted in relative units with v(4T)/J = 1 for H\\a, and 
v(6T)/J = 1 for H\\c. 



for H\\c. Comparing the fitting result to Eq. 3.6, 3.16, we obtain 



g » - g aa D/2J = 0.04 



(3.19) 



g s + g cc D/2J = 0.12, 



where g aa = 2.28 and g cc =2.12 are given by Eq. 3.1. The results, g s « 0.08 and D/J 
~ 0.04, are quite plausible; for comparison, D/J in Copper Benzoate and g s ~ 0.02 
and D/J w 0.1 [50]. 

We conclude that the quantum Sine-Gordon model explain the temperature de- 
pendence of specific heat and field induced gap quite well. We were able to determine 
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the field dependence of spin wave velocity, and estimate the values for DM vector, D, 
and staggered g tensor, g s . However, the gap values in the high field do not follow the- 

2 

oretical prediction A m H~z that well. This might be a result of a 3D ordering phase 
transition in CDC, which is not observed in Copper Benzoate. which is not observed 
either in Copper Benzoate[37], in [PM • Cu(N0 3 ) 2 ■ (H 2 0) 2 ]„[96], or in Yb 4 As 3 [97]. 

3.5.2 Neutron diffraction 

To resolve the magnetic ordered structure of CDC, we first determine what are 
the magnetic structures compatible with the symmetry of the crystal. The symmetry 
analysis is based on the assumption that any magnetic structure adopted through 
a second order phase transition can be expanded in the basis functions for a single 
irreducible representation of the magnetic space group. The space group of CDC is 
G=Pnma (D 2 h 16 ) [98]. There are four magnetic Copper ions in the unit cell occupying 
the positions 

4(c) : l(x, 1/4, z),2(x + 1/2, 1/4, 1/2 - z), 3(1/2 - x, 3/4, z + 1/2), 4(1 - x, 3/4, 1 - z), 

(3.20) 

where x= 0.309 and £=0.3614 [98]. Each of the copper spins have three components. 
The 12 components form a representation for the magnetic structure and group theory 
can be used to reduce the representation to its irreducible representations. Because 
magnetic Bragg peaks coincide with forbidden nuclear peaks, the wave vector associ- 
ated with the magnetic structure should be K=0. Using the standard group theory 
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analysis[104], the magnetic representation can be decomposed into irreducible repre- 
sentations as follows 

d M k = n + 2r 2 + 2t 3 + r 4 + r 5 + 2r 6 + 2r 7 + r 8 . 

The basis functions for the irreducible representations have been calculated as shown 
in Table 3.4. The actual spin directions can be combinations of the basis func- 
tions within a single irreducible representation. The magnetic structure is deter- 
mined by comparing measured and calculated integrated magnetic Bragg peak in- 
tensities of all the basis functions listed in Tab 3.4. The basis function which pro- 
vides the best account of measured intensities, i.e., the smallest residual coefficient 
R = Eq l /o6s (Q) - /ca '(Q)l/Eq /o6s (Q) describes the magnetic structure of CDC. 
To calcuate R, magnetic Bragg peak intensities are brought to the absolute scale by 
comparison with nuclear Bragg peak intensities measured in the same configureation. 
Table 3.6 lists the integrated intensities of nuclear Bragg peaks observed in the ex- 
periment and calculated from known crystal structure. They agree with each other 
reasonably well. Table 3.5 gives a list of the calculated residual coefficients for all 
the irreducible representations. At zero field and T = 0.25K, the smallest residual 
coefficient R=0.25 is given by representation: s'm9 r 2 + cos 9 t' 2 with 6 = 17° ± 10°. 
The ordered moment determined from this structure was \M\ = 0.741/_i B per occupied 
copper site. Fig. 3.14 (a) shows a schematic view of the spin arrangement of the the 
ordered phase with spin vectors in ac plane. The magnetic structure in the spin flop 
phase is also determined. The magnetic structure given by representation T4 provides 
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Table 3.4: The basis functions for irreducible representations of space group Pnma 
(DID w ith magnetic wave vector K=0 and for atoms on the 4c site. 



the best account of the measured intensity at H = 0.5T and T = 0.25K data with 
R=0.188 and \M\ = 0.726fi B . The ordered spin vectors are collinear and along the b 
axis, as shown in Fig. 3.14 (b). The calculated intensities of magnetic Bragg peaks 
from the spin configuration described above are listed in Table 3.7 with the measured 
intensities, they are in good agreement with each other. 



3.5.3 Conclusion 

In summary, we have presented that the heat capacity and neutron diffraction 
measurements of S=l/2 antiferromagnetic chain system CDC. A 3D ordered phase 
has been observed and the ordering temperature as a function of the magnitude of 
the field and its orientation has been determined. There is a field induced gap due 
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irreducible 


TT 

n 


= 0i 


ti = 


0.5i 


representation 


Rmin 


m(fjL B ) 


Rmin 




Tl 


0.905 


0.464 


0.915 


0.375 




0.250 


0.741 


0.681 


0.643 


r 3 


0.670 


0.755 


0.83 


0.470 




0.358 


0.77 


0.188 


0.726 


T 5 


1.00 




1.000 




T 6 


0.725 


0.786 


0.957 


0.258 


T7 


0.928 


0.342 


0.982 


0.170 


Tg 


1.00 




1.00 





Table 3.5: Residual coefficients i? = E r \I ohs {r) - I cal {t)\ j ^ r 7 o6s (r) were calculated 
for all the magnetic structure described by the basis functions. The magnetic struc- 
ture of CDC is described by the basis function which provides the best account of 
measured intensities, i.e., the smallest R value. 
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Table 3.6: The measured and calculated integrated intensities of nuclear Bragg peaks 
for CDC. They agree with each other reasonably well. 
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Figure 3.14: A schematic view of magnetic ordered structure of CDC showing only 
the copper atomic sites. The atoms are labeled as in Eq. 3.20. Vectors represent 
the ordered spin structures as determined by elastic neutron scattering, (a) Ordered 
structure in zero field, spin vectors are in the ac plane with = 17° ±10°. (b) Ordered 
structure in the spin flop phase, spin vectors are along b axis. 
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Table 3.7: The measured and calculated Q-integrated intensities of magnetic Bragg 
peaks of CDC. Magnetic Bragg intensities were obtained from the difference between 
the scattering intensity at 0.25 K (well below Neel temperature) and 2 K (above Neel 
temperature). Im 1 { t ) was calculated using the spin configuration described in the 
text. 



to an effective staggered field and the gap scales approximately as a 2/3 power of 
the applied magnetic field. This is in good agreement with the quantum sine-Gordon 
Model. We also determined the magnetic structure and magnetic moment of the 
ordered phase. Future experiment will focus on the inelastic neutron scattering under 
an applied magnetic field to map soliton and breather contributions to the dynamic 
spin correlation function. 



61 



Chapter 4 

S=l Heisenberg antiferromagnetic 
chain: NDMAP 

4.1 Introduction 

Analogous to the S=l/2 chain in a field, the spectrum of an isotropic S=l chain is 
gapless for H > H c and dynamic spin correlations are expected to have a low energy 
incommensurate mode, however, it is challenging to study magnetized spin-1 chains 
experimentally. This is because most of existing experimental systems can only be 
magnetized by a very high magnetic field due to the Haldane gap. Moreover, in 
some systems, local symmetry axis alternate along spin chains leading to a staggered 
field which induces a gap in the spectrum and possibly disfavors incommensurate 
soft modes. However, the recent discovery of NDMAP (Ni(C 5 Hi4N2)2N 3 (PF 6 ))[105], 
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a relatively easily crystallized Haldane-gap compound with a critical field of only 
around 4 T[74], and no staggering of g-tensors within spin chains, promises to make 
the high-field phase accessible to inelastic neutron scattering studies. 

4.2 NDMAP 

The crystal structure of NDMAP is orthorhombic (space group Pnmn, a = 
18.046 A, b = 8.705 A, and c = 6.139 A) as shown in Fig. 4.1. The AF spin chains 
run along the c axis and are composed of octahedrally-coordinated 5 = 1 Ni 2+ ions 
bridged by triplets of nitrogen atoms. These long nitrogen bridges account for a rela- 
tively small in-chain AF exchange constant J = 2.6 meV, as estimated from bulk x(T) 
measurements [74]. Specific heat studies revealed a transition to the high-field phase 
at He = 3.4 T, and H^ = 5.8 T, extrapolated to T — > 0, for a magnetic field applied 
parallel and perpendicular to the chain-axis, respectively [74, 106]. This anisotropy of 
critical field is attributed to single-ion easy-plane magnetic anisotropy of type DS 2 Z . 
The anisotropy constant was obtained from bulk susceptibility data: DjJ m 0.3. 
Using the numerical result A z « 0.41J + 2pD, A x w 0.41J - pD, p w 2/3[61, 83], 
one can estimate the Haldane g ctp energies! ^^|| ~ 

2.1 meV and f« 0.54 meV. 



63 




Figure 4.1: A schematic view of the antiferromagnetic spin chains in the NDMAP 
crystal structure. Only Ni, N and C atoms are shown. 



4.3 Spin Hamiltionian of NDMAP 
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4.3.1 Introduction 

While NDMAP appears to be an ideal model system for neutron scattering experi- 
ments in the high-field phase, additional characterization, particularly a measurement 
of exchange constants and and anisotropy parameters, is required before such a study 
can be carried out. Finding the 3D AF zone-center is especially important, since it is 
at this wave vector that static long-range correlations are expected to appear in the 
high-field phase. We performed zero-field inelastic neutron scattering experiments on 
deuterated NDMAP single-crystal to obtain this information. 

4.3.2 Experimental procedures 

Most inelastic neutron scattering data were collected on a 140 mg deuterated sin- 
gle crystal NDMAP sample. It was observed that the crystals tend to shatter when 
cooled to low temperature, and even more so when warmed back up. When wrapped 
in aluminum foil the crystals do not fall apart, but the width of the mosaic distri- 
bution increases with thermal cycling (Fig. 4.2). So we collected data on the single 
crystal sample that was taken through the cooling cycle no more than twice. The 
mosaic of the as-grown sample was roughly 0.3° and increased to 1.5° and 3° after 
the first and second cooling, respectively. High-resolution inelastic neutron scattering 
measurements were performed on the SPINS triple-axis spectrometer installed on the 
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Figure 4.2: Rocking scans through a nuclear Bragg peak (2, 0, 0). At room tempera- 
ture, the full width at half maximum is 0.40(1)°. After first cooling down to 0.25 K, 
the full width at half maximum increases to 0.99(2)°. 
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cold neutron source at the National Institute of Standards and Technology Center for 
Neutron Research (NCNR). Pyrolytic graphite crystals set for their (002) reflection 
were used as monochromator and analyzer. Beam divergences were approximately 
40' — 80' — 80' — 240' through the instrument, with a cooled Be filter between the 
monochromator and sample. The measurements were done with a fixed final neutron 
energy Ef = 2.8 meV. The crystal was mounted with either the a or b crystallo- 
graphic axis vertical, making (0,k,l) and (h, 0, I) reciprocal-space planes accessible 
to measurements. The sample was cooled to T = 1.4 K in an "ILL-Orange" cryostat 
with a 70 mm diameter sample well. 

4.3.3 Experimental results and data analysis 

Most of the data were collected in constant-Q scans at the ID AF zone-center / = 
0.5 in the energy transfer range 0-1 meV. In these measurements of the 3D dispersion 
relation, wave vector resolution was particularly important, so a flat analyzer was 
used. Typical scans for different momentum transfers perpendicular to the c-axis 
are shown in Fig. 4.3. At all wave vectors a well-defined peak corresponding to 
the lower-energy Haldane excitation doublet is clearly seen around 0.5 meV energy 
transfer. To observe the higher-energy Haldane-gap excitation and obtain an accurate 
measurement of the anisotropy constant we performed a constant-Q scan in the range 
0-2.4 meV (Fig. 4.4). To maximize intensity we used a horizontally-focusing analyzer 
pointing c* towards the analyzer to maintain wave vector resolution along the chain. 
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Figure 4.3: Typical constant-Q scans measured at T = 1.2 K in an NDMAP deuter- 
ated single-crystal. Solid lines are based on global fits to the data, as described in 
the text. The shaded areas represent partial contributions of the o-axis and 6-axis 
polarized Haldane gap excitations to the intensity. The dashed line indicates the 
background level. The magnetic excitations are resolution-limited and peak shapes 
are entirely defined by resolution effects. 
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The scattering angle was varied so the projection of wave vector transfer along the 
chain was 0.5c* throughout the scan. In addition to the peak seen in the low-energy 
scans, a weaker feature is observed at hu ~ 2 meV that can be attributed to the 



c-axis-polarized Haldane gap mode. 



h (r. 1. u.) 



4.0 
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h(£> (meV) 

Figure 4.4: A constant-Q|| scan measured in NDMAP using a horizontally-focusing 
analyzer configuration. Top axis shows the momentum transfer perpendicular to the 
chain-axis. Lines are as in Fig. 4.3. 



The steep dispersion along the chain direction is best measured in constant-E 
scans well above the gap energy. Given the small sample size, intensity, which for 
antiferromagnetic spin waves scales as 1/co, becomes a serious limiting factor. To 
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overcome this problem, the measurements were performed in a coarser resolution 
setup using a thermal neutron beam at the BT-2 triple-axis spectrometer at NCNR. 
14.7 meV fixed-incident energy neutrons were used with a PG filter in front of the 
sample and 60' — 60' — 60' — open' collimations. Constant-energy scans were collected 
at 3 meV and 4 meV energy transfer near the ID AF zone-center I = 1.5 (Fig. 4.5). 
The two spin wave peaks are clearly resolved in both scans. The higher-/ peak appears 
sharper than the lower-/ feature due to focusing effects. Note that only transverse 
excitations contribute to scattering in this geometry. 

The measured constant-Q scans were analyzed using a parameterized model cross 
section, numerically convoluted with the Cooper-Nathans triple-axis spectrometer 
resolution function [107]. Near the ID AF zone-center the single mode approximation 
(SMA) for the dynamic structure factor of isolated Haldane spin chain is known 
to work extremely well[65, 71]. For each channel of spin polarization the dynamic 
structure factor 5^ aQ ^(q, uj) can be written as: 

SM&u) oc |/(,)| 2 f 1 - 1 " T (QC) ^5(hu> - ^ Q , q ), (4.1) 

(^ Q , q ) 2 = A 2 + u 2 sin 2 (qc). (4.2) 

Here v is the spin wave velocity, given by v ~ 2.49J[108]. The dimensionless con- 
stant Z defines the static staggered susceptibility of a Haldane spin chain: x-k = 
Zv/ A..(Reie we have adopted the notation used in Ref. [109]) In the above expres- 
sion we have included the magnetic form factor f(q) for Ni 2+ and the polarization 
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Figure 4.5: Constant-energy scans near the ID AF zone-center I = 1.5 (raw data) , 
showing the dispersion of Haldane-gap excitations in NDMAP along the chains-axis. 
Solid lines are a fit to the data as described in the text. 
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factor ^1 — ^f^J . The structure factor for weakly-coupled chains can be calculated 
in the Random Phase Approximation (RPA)[79]. The expression for S( a °)(q,aj) does 
not change explicitly, but the excitations acquire dispersion perpendicular to the chain 
axis: 

(hcu^) 2 = A 2 a + v 2 sin 2 (qc) + ZvJ'(q). (4.3) 

In this formula J'(q) is the Fourier transform of inter-chain exchange interactions, 
which we assume to be isotropic (of Heisenberg type). According to numerical cal- 
culations Z 1.26 [108]. The form of J' can be guessed by looking at the crystal 
structure. The smallest inter-chain Ni-Ni distance (8.705A) is along the b crystal- 
lographic axis. We shall denote the corresponding exchange constant as J y . The 
next-smallest inter-chain Ni-Ni distance (10.478A) is along the (0.5,0.5,0.5) direction. 
This interaction, however, is frustrated by in-chain AF interactions (any site in one 
chain couples to two consecutive sites in another chain), and is thus irrelevant, within 
the RPA, at momentum transfers qc ~ tt. Finally, the third-smallest inter-chain 
distance (18.046A) is along a. The corresponding coupling constant J x is expected to 
be very small, due to the large atomic spacing. The Fourier transform of inter-chain 
interactions can thus be written as: 

J'(q) = 2 J x cos(qa) + 2J y cos(qb). (4.4) 

The relevant parameters of the model are thus the three gap energies A x , A y and 
A z , three exchange parameter J, J x and J y and an overall intensity prefactor. For 
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constant-Q scans three additional parameters were used to describe the background: 
an energy-independent component, and the intensity and width for a Lorentzian pro- 
file centered at hco = to account for incompletely resolved elastic scattering. The 
background was assumed to be flat in the constant-^ scans measured in the thermal- 
neutron setup. 

The difference between and A y , usually a result of in-plane magnetic anisotropy 
of type E(S% — Sy), is too small to observe two and separate corresponding peaks 
in any single constant-Q scan. Distinct polarization factors for the two branches 
however allow us to extract both parameters in a global fit to the data measured 
at different momentum transfers perpendicular to the chain-axis. As a first step, 
we simultaneously analyzed the data collected in (0,k,l) reciprocal-space plane (11 
scans with a total of 382 data points). A very good global fit using Eqs. 4.1 and 4.3 is 
obtained with a residual x 2 = 1-6- The solid lines in Fig. 4.3a-c were calculated from 
the globally optimized parameter set. The shaded areas represent the contribution 
of each mode, and the dashed lines show the background level. Figure 4.6 (right) 
shows the obtained 6-axis dispersion relation (solid line), that has a minimum at 
k = 0.5. Symbols in this figure represent the excitation energies obtained in fits 
to individual scans. A similar global fit (4 scans, 149 data points, \ 2 = 1-7) was 
applied to all constant-Q scans measured in the (h, 0, /) scattering plane (solid lines in 
Fig. 4.3d and e). In this case the excitation energies at the zone-boundary (0, 0, 0.5) 
were fixed to the values obtained from the (0, k, Z)-plane global fit. Only the a- 
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axis exchange constant J x was refined. Dispersion of magnetic excitations along this 
direction is barely detectable, J x = being within the experimental error bar. The 
resulting dispersion relation is nontheless plotted in a solid line in Fig. 4.6 (left). 
The final values of parameters determined in the least-squares analysis are as follows: 
A x = 0.42(0.03) meV, A y = 0.52(0.06) meV, J y = 1.8(0.4) • 10" 3 meV, and J x = 
3.5(3.0) • 10~ 4 meV. From the inelastic measurements we can thus only conclude that 
the global minimum of the 3D dispersion (3D AF zone-center) is located somewhere 
on the (h, 0.5, 0.5) reciprocal-space rod, with only a hint that (0, 0.5, 0.5) is the actual 
critical wave vector. 

To verify that the splitting of the excitation doublet is indeed a relevant parameter, 
we have also analyzed the measured constant-Q scans assuming A x = A y . To obtain 
good fits to individual scans using this model, we were forced to allow the magnetic 
excitations to have a non-zero intrinsic energy width. However, even including this 
additional parameter did not yield to a good global fit to all measured constant-Q 
scans. The excitation energy determined from fits to scans collected in the (0, k, I) 
plane at large values of momentum transfer perpendicular to the chain axis were 
found to be smaller than at equivalent reciprocal-space points at smaller k. This 
behavior is consistent with having two separate resolution-limited modes at slightly 
different energies, the intensity of the higher-energy mode being suppressed by the 
polarization factor at large k. The data is inconsistent with the A x = A y model. 

The energy gap for c-axis polarized excitations was determined by fitting the 



74 




Figure 4.6: Transverse dispersion of lower-energy Haldane-gap excitations in NDMAP, 
determined in global fits to constant-Q scans (solid lines), as described in the text. 
The data points represent excitation energies obtained in fits to individual scans. 
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model cross section to the wide-range constant-Q scan measured with the horizontally- 
focusing configuration (Fig. 4.4, solid line). In this procedure the values of and 
A y were fixed. The analysis yields A z = 1.86(0.1) meV. 

The intra-chain exchange constant was determined by fitting the SMA cross- 
section to the measured const-i£ scans. The gap energies were fixed at the values 
quoted above and only J, a flat background and an intensity factor were refined to 
best-fit the data (solid line in Fig. 4.5. The obtained value J = 2.85(0.05) meV is in 
a good agreement with the result of Ref. [74]. The anisotropy constant D is straight- 
forward to estimate from the measured gap energies:[61, 83] D m |(2A 2 — A x — A y ) = 
0.70 meV, and Dj J = 0.25, which is also consistent with the bulk susceptibility result 
of Ref. [74]. The relative strength of inter-chain interactions are J y j J ~ 5.6- 10~ 4 and 
J x /J « 1.0 • 10~ 4 . These ratios are very similar to those found in NENP[71, 72, 73]. 

4.3.4 Discussion 

While the gap energies and the spin wave velocity measured with inelastic neutron 
scattering in this work are in excellent agreement with the results of susceptibility 
and high-field magnetization and specific heat studies [74], they appear to be a dis- 
crepancy with the ESR results of Ref. [106]. Sharp ESR absorption lines in NDMAP 
occur at 0.05 meV, 0.42 meV and 0.48 meV, extrapolated to zero field. For a uniform- 
chain Haldane antiferromagnet, low-energy ESR lines are associated with transitions 
between components of the Haldane triplet, that become thermally populated at 



76 

non-zero temperature. The energy of the 0.05 meV ESR line could correspond to 
a transition within the lower-energy doublet of transverse excitations (0.1 meV. ac- 
cording to our neutron data). In contrast, the other two resonance energies are much 
smaller than the doublet-singlet splitting observed in this work, and that can be ex- 
pected from the known relations between D, J, and A. In fact, the resonances are 
rather close to the values of the transverse Haldane gap energies themselves. However, 
in the absence of a staggering in the spin chains, ESR transitions from the ground 
state to single-magnon states should be prohibited by symmetry. Further studies will 
be required to fully understand the ESR spectra. 

The small value of critical fields in NDMAP, compared to those in NENP, are 
a result of smaller in-chain exchange interactions and somewhat larger easy-plane 
magnetic anisotropy. The lowest-energy excitation in NDMAP is polarized along the 
crystallograhic a-axis. The 3D magnetic zone-center is located on the (h, 0.5,0.5) 
reciprocal-space rod. 

4.4 Field-induced three and two-dimensional 
freezing 

4.4.1 Introduction 

The magnetization process of a S = 1 ID AF is now rather well understood 
theoretically[85, 110, 111, 112, 113, 114]. For the isotropic case the problem was 
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shown to be equivalent to Bose condensation in one dimension [85, 115]. Many of the 
theoretical results were confirmed in experimental studies of real quasi-lD compounds. 
For a long time however, the actual phase transition remained inaccessible for experi- 
mental investigation. In some materials (such as Y 2 BaNiOs, for example[63, 64, 65]), 
the value of H c is prohibitively high. In other compounds {e.g.. NENP [71, 72, 73]), 
the transition does not occur due to the staggered g-tensor, and is instead replaced by 
a broad cross-over phenomenon[84, 116, 117, 118]. A true phase transition at H c has 
been experimentally observed only recently in NDMAP [74, 106, 119] and NDMAZ 
[120]. Specific heat, magnetization and ESR studies have provided a comprehensive 
picture of the H — T phase diagram [74, 106, 119], a refined version of which shown 
in Fig. 4.7. Nevertheless, to date, the nature of the high-field phase remained unde- 
termined, and no direct information about static magnetic structure in the high field 
phase has been obtained experimentally. Here we present the results of high-field 
magnetic neutron diffraction studies of NDMAP. We find that the high-field phase 
is characterized by a commensurate ordering of spin components perpendicular to 
the field direction. Surprisingly, the nature of the high-field phase depends on the 
direction of applied field: depending on the experimental geometry, for H > H c the 
system is either in a Neel state with true 3D LRO, or a spin glass phase with highly 
anisotropic spin correlations. 
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Figure 4.7: Field-temperature phase diagram of NDMAP. A quantum-disordered 
(QD) phase with a gap is seen in low fields. The high-field phase is characterized 
by static long-range order (3D-LRO) or quasi-2D short-range order (2D-SRO). Open 
symbols from specific heat measurements, as in Ref. [106]. Solid symbols from neutron 
diffraction. Lines are guides to the eye. 
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4.4.2 Experimental procedure 

The neutron diffraction experiments were performed at the NG5-SPINS cold- 
neutron 3-axis spectrometer at the National Institute of Standards Center for Neutron 
Research. A 150 mg sample was mounted with either the (0, k, I), (h, h, I) or (h, k, k) 
reciprocal-space planes in the scattering plane of the spectrometer. Most measure- 
ments were performed using (guide) — 80' — 80' — (open) or (guide) — 40' — 40' — (open) 
collimations, with a Be-filter positioned in front of the sample to eliminate higher- 
order beam contamination, and Ef = 5 meV or Ef = 3 meV fixed final energy neu- 
trons. Sample environment in all cases was a pumped 3 He cryogenic insert and a 9 T 
or 7 T superconducting magnet, with the field applied vertically, i.e., perpendicular 
to the scattering plane. 

4.4.3 Experimental results and data analysis 

Three dimensional commensurate AF LRO was observed in a crystal mounted in 
the (h,k,k) zone, with the magnetic field applied along the (0, —1, 1) direction. At 
T = 0.25 K the transition occurs at H c = 4.8 T, and manifests itself in the appearance 
of new Bragg reflections, at (h, ^y^, ^y^) (h-even integer; n- integer) reciprocal- 
space positions. No additional peaks were found at odd-h reciprocal-space points. 
Figures 4.8(a) and 4.8(b) show typical scans through the (0, 0.5, 0.5) peak measured 
at H = 9 T applied field (symbols). The peaks were found to be resolution-limited 
along both the (1, 0, 0) and (0, 1, 1) directions. (A very small intrinsic width of the 
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Figure 4.8: Typical elastic scans (open symbols), showing magnetic Bragg scat- 
tering in NDMAP at T = 0.25 K, H = 9 T|| [0,-1,1] (a,b) and T = 0.375 K, 
H = 7 T||[1,0,0] (c-e). The solid lines and symbols represent experimental wave 
vector resolution. Dashed lines as described in the text. 
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peak Ri 0.002 ± 0.0004 A -1 is extracted by fitting the scan along (0. 1,1) direction to a 
Lorentzian profile convoluted with the resolution function.) Experimental resolution 
(Fig. 4.8, solid lines) was either calculated or directly measured in our experiment. 
Geometrical limitations imposed by the use of the bulky 9 T magnet prevented us 
from measuring the width of the peak in the vertical (0, 1, —1) direction. However, 
the substantial values of the b and c axis coupling constants, as well as the resolution- 
limited peak widths along a, where magnetic interactions are weakest, suggest that 
the observed reflections indeed have zero intrinsic width in all three directions, and 
represent true long-range spin correlations in 3 dimensions. This is also consistent 
with the measured absolute values of magnetic intensities, as discussed below. To 
determine the spin arrangement in the high-field ordered phase, intensities of the (0, 
0.5, 0.5) (2, 0.5, 0.5) (4, 0.5, 0.5) (2, 1.5, 1.5) and (4, 1.5, 1.5) magnetic Bragg peaks 
in the (h,k,k) plane were measured in rocking scans at T = 0.25 and H = 9 T. 
The magnetic intensities were normalized to an absolute scale by comparing them 
to nuclear intensities measured in the same configuration at room temperature (the 
exact low-temperature structure of NDMAP has not been determined to date). The 
data were analyzed using a simple collinear model. As the energy scale of in-chain 
AF interactions is much larger than that defined by the applied magnetic field ( J 
gfiBH), the ordered moment is expected to be perpendicular to the field direction. 
The (a, 6)-easy plane anisotropy ensures that it should also be perpendicular to the 
chain-axis. For a field applied along [0,-1,1] the spins are thus along (1,0,0), as 
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confirmed by recent ESR experiments[106]. This allows only one magnetic structure 
that would be consistent with the observed selection rule for magnetic Bragg peaks. 
Varying only the sublattice magnetization m, we obtained a good fit to the data 
(x 2 = 1.6) with m = 1.13(5) /ie- This value is to be compared to m = 2 fj, B 
for a fully saturated S = 1 antiferromagnet. The measured field dependence of the 
(0, 0.5, 0.5) magnetic Bragg peak intensity at T = 0.25 K is shown in Fig. 4.9(a) (solid 
circles), and is consistent with our expectations for an easy-plane Haldane system at 
T — > 0. The data were analyzed assuming a power-law field dependencies of the 
Bragg intensity and ordered moment: 1(H) oc (H — H c ) 2 ^ . A fit of this equation 
to the data (Fig. 4.9(a), solid line) yields H c = 4.81(1) T and j3 = 0.207(5). Since 
the transition at H c occurs through a softening of the lowest-energy Haldane gap 
excitation [110], H c is determined by the Haldane gap energies. It is straightforward 
to show that for a magnetic field applied at angle a to the anisotropy axis c is given 
by: HbH c = A z A±/ ^fA| cos 2 a + g 2 L A z A± sin 2 a. In this formula g z = 2.1 and 
g± = 2.17 are components of the Ni 2+ gyromagnetic tensor along and perpendicular 
to the anisotropy axis, respectively [74]. For H||[0, —1, 1], a = 54.8° and this equation 
gives H c = 5.4 T, in reasonable agreement with our experimental result. In the case of 
broken axial symmetry, for a single chain in a magnetic field at T = 0, the transition 
is expected to fall in the 2D Ising universality class, with magnetic field taking the 
role of effective temperature[85]. The order parameter critical exponent for this model 
is (3 = 0.125. For NDMAP however, inter-chain interactions along the b axis can be 
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considered substantial, and 3D Ising behavior, with (3 ~ 0.31, may be expected. The 
measured critical exponent falls in between these two values indicating a dimensional 
crossover regime. 

The temperature dependence of the (0,0.5,0.5) magnetic Bragg intensity mea- 
sured at H = 9 T is shown in Fig. 4.9(b) (solid circles) . An analysis of the data 
assuming a power law (solid line) yields T c = 1.7(1) K. The critical exponent was 
found to be indistinguishable from the mean field value /3 = 0.5, indicating that at 
T k 1.5 K the true critical region is too narrow to be investigated in the present 
experiment. A striking result of this work is that for a magnetic field applied per- 
pendicular to the chain-axis, the high-field phase is no longer a 3D-ordered state, 
but has predominantly 2-dimensional spin correlations. For the crystal mounted in 
the (0,k,l) zone, at T = 0.38 K and H|| [1,0,0], magnetic scattering was detected 
at (0, ^y^, 2m 2 +1 ) (n, m-integer) reciprocal-space points above H c = 6 T. k- and l- 
scans through the (0, —0.5, 1.5) reflections contain well-defined peaks, as shown in 
Figs. 4.8(d) and 4.8(e) (symbols). Scans perpendicular to the (b,c) plane however, 
reveal that the scattering is concentrated in Bragg rods along the (1,0.0) direction, 
rather than Bragg peaks, as in the H||[0, —1, 1] case. Figure. 4.8(c) shows an /i-scan 
through the (0, —0.5, 1.5) position, that could only be performed in a limited range 
due to geometrical constraints. The measured intensity is independent of h. Note 
that the measured vertical resolution (Fig. 4.8(c), solid line), was sufficient to observe 
a well-defined Bragg peak. The rod nature of magnetic scattering implies that static 
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spin correlations along the a axis are absent in the system, and that magnetic or- 
dering occurs only within individual (b, c) planes. This is consistent with the a-axis 
being the direction of weakest magnetic interactions. Similar behavior was observed 
for a magnetic field applied along the [1,-1,0] direction, where magnetic scattering 
was observed above H c = 5.8(2) T, at T = 0.25 K, and also found to be concentrated 
in Bragg rods parallel to the a axis. In this geometry, the rods cross the (h, h, I) scat- 
tering plane at (^^, 2m 2 +1 ) {n, m- integer) positions. The measured field- and 
temperature dependencies of magnetic diffraction intensity at q = (0.5,0.5,0.5) are 
shown in Fig. 4.9 (open symbols). A well-defined transition is observed in agreement 
with specific heat measurements (Fig. 4.7, diamonds). 

The apparent Bragg rod intensity for H||[l, 0, 0] and H||[l, —1. 0] is substantially 
smaller than the intensity of magnetic Bragg peaks seen for H||[0, 1, —1] (Fig. 4.8). 
This is due to that fact that in the 2D-ordered phase magnetic intensity is spread out 
through an entire Brillouin zone along the a axis, rather than being concentrated at 
h = 0. However, in the 2D phase, within each (6, c) plane, the actual ordered moment 
is similar to that in the 3D ordered state. The 2D ordered moment for H||[1,0,0] 
and H||[l,— 1,0] was deduced from an analysis of measured Bragg rod intensities, 
paying close attention to resolution effects associated with vertical angular acceptance 
of the spectrometer. A series of rocking scans across the rods were collected for 
H = 7 T||[1,0,0], T = 0.38 K and H = 9 T||[l,-1,0], T = 0.25 K, respectively. 
The data were normalized by the measured nuclear Bragg peak intensities. The 2- 
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H(T) T(K) 

Figure 4.9: Measured field (a) and temperature (b) dependence of the ordered moment 
in NDMAP for a magnetic field applied in the (b, c) plane. Solid lines are power- 
law fits to the data, as described in the text. Inset: A schematic view of the spin 
arrangement in the high-field phase, for the case of H||[0, -1, 1]. 
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dimensional spin structure was assumed to be collinear, with spins perpendicular to 
both the chain axis and H in each case. Nearest-neighbor spins were assumed to 
be aligned antiparallel along both the b and c axes. This simple model was found 
to agree very well with the available data. The ordered moment was determined 
to be m = 0.66(4) fi B and m = 1.2(1) fx B for # = 7 T|| [1,0.0], T = 0.38 K and 
H = 9 T||[l,— 1,0], T = 0.25 K, respectively. The latter value agrees very well 
with the 3D ordered moment measured at H||[0, —1, 1] at the same value of field and 
temperature. 

A careful analysis of the scans across the Bragg rods for H| | [1, 0, 0] reveals a finite 
correlation length even within the quasi-2D ordered (b, c) planes. Along the direction 
of strongest coupling, i.e., along the chains, the magnetic peaks are resolution-limited. 
This can be seen from Fig. 4.8(e), where the solid line represents the experimental 
Z-resolution. However, along the b axis, where the magnetic interactions are weaker, 
the scattering rods are visibly broader than experimental resolution (Fig. 4.8(d), 
solid line). Fitting the &-scan to a Lorentzian profile convoluted with the resolution 
function (Fig. 4.8(d), dashed line) allows us to extract the intrinsic width of the rod 
that corresponds to a real-space correlation length £(,=100(20) 106. The high-field 
phase in this case is thus characterized as a quasi-two dimensional ordered state with 
almost-perfect correlations along the chains, medium-range correlations along the b 
axis and only very short range correlations along the a axis. 
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4.4.4 Discussion 

The type of "spin freezing" in NDMAP is very similar to that recently found in 
the 5=1/2 quasi-lD AF SrCu02 [19]. For a two-dimensional frustrated Heisenberg 
system, the ground state is disordered and gapless for half-integer and odd-integer 
spin due to "instanton" topological defects [122]. If magnetic coupling along the 
third direction is sufficiently weak, as it maybe in NDMAP, it will not be possible 
to restore LRO in the system due to pinning of such defects within each plane. The 
similar behavior of two very different systems, NDMAP and SrCu0 2 , suggests that 
anisotropic spin freezing may be a rather general feature of quasi-low dimensional 
antiferromagnets. What makes the case of NDMAP special is that both true 3D 
LRO and spin freezing can occur in the same sample, depending on the direction of 
applied field. This may result from the anisotropic nature of magnetic interactions 
along the a axis, that, due to the very long lattice spacing, are expected to be largely 
dipolar. The strength and even sign of such coupling depends on the orientation of the 
ordered moment, which, in turn, is defined by the direction of applied field. Changing 
the field orientation for NDMAP is thus a way to tune those magnetic interactions 
ultimately responsible for 3D LRO or spin-freezing behavior. 

For the basic physics of Haldane spin chains, a very important experimental result 
is the observation of a commensurate ordering of transverse (relative to H) spin 
components for fields above H c . Such staggered LRO in a quasi-lD material is a direct 
consequence of a divergence in the transverse spin correlation function for an isolated 
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chain in an external field at q = n. Theory predicts that above H c a divergence also 
exists in the longitudinal correlator, but at a field-dependent incommensurate wave 
vector [85, 113]. No incommensurate magnetic peaks were found in our experiments 
on NDMAP. Of course, a divergent susceptibility in a single chain does not guarantee 
the appearance of LRO at the same wave vector in a 3D material. It is also conceivable 
that incommensurate elastic scattering in NDMAP appears at a different wave vector 
transfer perpendicular to the chain axis, due to anisotropy of inter-chain interactions, 
and thus has escaped detection in the present study. Further work will be required 
to fully resolve this issue. 

In summary, our results provide direct experimental evidence of field-induced com- 
mensurate AF order in a Haldane-gap system, and are in good agreement with theo- 
retical expectations. Specifics of inter-chain interactions in NDMAP result in a rich 
phase diagram with the high-field phase being either a true 3D ordered state or a 
highly anisotropic short-range ordered state. 

4.5 Quasi-elastic scattering in the high field phase 

4.5.1 Introduction 

Haldane-gap systems exhibit numerous unusual properties, and particularly inter- 
esting predictions were made for their behavior in high magnetic fields [85, 110, 114, 
123, 126, 127, 128, 129]. In an external field the gap excitations, which are a S = 1 
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triplet, are subject to Zeeman splitting. The gap A for one of the three branches 
decreases with field (Fig. 4.10(a)) and closes at some critical field H c ~ A/g/_i B 
[84, 85, 110. 126. 127]. Near H c this problem is equivalent to Bose condensation 
[85, 110, 127]. However, at H > H c the formation of a condensate is prevented by 
strong interactions between magnons. As a result, the high-field phase is rather un- 
usual, with power-law spin correlations, and is an example of a Luttinger spin liquid 
[128, 129]. The basic physics at H > H c can be understood within a simple model. 
The magnon interactions are assumed to be a hard-core repulsion, in which case the 
excitations behave as free fermions. [85, 110, 124, 125, 127]. At H > H c the ground 
state is a Fermi sea of excitations with field-dependent Fermi density and wave vector 
qp, as illustrated in Fig. 4.10(b). The spectrum is gapless, and the dynamic spin 
correlations are incommensurate. The effect is related to field-induced incommensu- 
rability in gapless spin chains, that is described using a different fermion mapping 
[12, 37, 130]. Ni(C 5 H 14 N2)2N3(PF 6 ) (NDMAP) has a critical field of only 5 T[74]. 
For H > H c residual interactions between the S = 1 Ni 2+ -chains in this compound 
lead to commensurate long-range AF ordering at low temperatures as discussed in 
Section 4.4 [74, 76]. However, as reported in Section 4.3, the strength of effective 
inter-chain coupling is about 0.2 meV, which is small compared to the in-chain ex- 
change constant J m 2.8 meV [77, 75]. At hco > 0.2 meV one can thus directly study 
the exotic dynamics of the ID model. In the present work we report inelastic neutron 
scattering studies of low-energy excitations in NDMAP, for magnetic field below and 



90 



3 



\ H<H 


/ 

c / 

A>0 


(a) 





H>H 

c 




(b) 


/a<o 



71 



71 



Figure 4.10: Dispersion in the lowest-energy Haldane-gap mode below (a) and above 
(b) the critical field H c (blue). In the free-fermion model, excitations with a negative 
energy at H > H c form a Fermi sea (cyan) with a characteristic Fermi vector q F 
Lowest energy excited states are viewed as "particles" for q it > q F or "holes" for 
\q\\ — 7r | < q F (red line). 



well above H c , at temperatures above the 3D-ordered phase. 



4.5.2 Experimental procedure 

Experiments were carried out on the SPINS triple-axis spectrometer installed 
at the cold neutron facility at the National Institute of Standards and Technology 
Center for Neutron Research. 23 deuterated NDMAP single crystals were co-aligned 
to produce a sample of total mass 1.2 g (Fig. 4.11) and a mosaic of about 7° (Fig. 
4.12. The sample was mounted in a cryomagnet and measurements were performed 
at T = 2.5 K and magnetic fields H = 0, 2, 4, 5.5, 7 and 9 T. The field was applied 
along the b axis. The spectrometer employed a vertically focusing pyrolytic graphite 
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(PG) (002) monochromator and a 23 cmxl5 cm flat multi-blade PG(002) analyzer 
positioned 91.4 cm from the sample. Scattered neutrons were registered by a position- 
sensitive detector (PSD). Measurements were done with fixed-final neutron energy for 
the central analyzer blade of 3.1 meV. The collimation setup was ( 58 Ni-guide)-(open)- 
(open)-80' (radial), and a Be filter was used either before or after the sample. In order 
to minimize the contributions of the sample mosaic to wave vector (Q) resolution, all 
scans were performed such that momentum transfer for the center pixel of the PSD 
was along the (00/) direction in reciprocal space. Neutrons registered by the detector 
correspond to scattering events in the sample with almost equal momentum transfers 
along the spin chains. In a single setting of the spectrometer the PSD covered an 
energy range of 1 meV and a range of wave vector transfer along the chain of 0.3 r. 1. 
u. Separate scans were corrected for the analyzer efficiency, and merged to produce 
a single data set. A linear background was subtracted from each const-i£ slice taken 
from the data set at each field. The data were smoothed using a Gaussian kernel with 
0.05 meV FWHM in energy and 0.01 r. 1. u. in /. 

4.5.3 Experimental results and data analysis 

Typical resulting data sets are visualized in Fig. 4.13. For the following discus- 
sion it is important to estimate the critical field H c . Long-range ordering, which is 
essentially a three-dimensional (3D) effect and occurs at H % C B m 12 T for T = 2.5 K, 
is of little interest in the present study. The relevant quantity is the field H c at which 
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Figure 4.11: (a) NDMAP crystals (b)Neutron scattering sample holder with 23 
NDMAP crystals inside. 
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Figure 4.12: Rocking scan through a nuclear Bragg peak (0, 0, 2) of the sample made 
of 23 coaligned crystals at room temperature. The width of the peak is about 7 °. 
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Figure 4.13: False color plots of inelastic intensity measured in NDMAP at T = 2.5 K 
for different values of magnetic field applied along the b axis. The extra intensity 
seen below q\\ = 0.4c* at low energies in (c) is most likely an artifact of imperfect 
background subtraction. Solid lines are as described in the text. 
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the gap for the lower-energy Haldane excitation vanishes at q = (0,0,0.5), where 
all the measurements were performed. H c can be estimated from the gap energies 
measured in zero field [77] and the known gyromagnetic ratio for the S=l Ni 2+ ions 
in NDMAP [74]. Using Eq. 2.14 in Ref.[83] we get H c « 7.3 T. Below we shall first 
review the results obtained for H < H c , and proceed to discuss the 9 T data, that 
clearly represent the high-field phase and contains the most important findings of this 
work. 

As established in Section 4.3, the Haldane triplet in NDMAP is split by single- 
ion anisotropy. The parabolic shape seen in Fig. 4.13(a) for H = corresponds to 
two lower-energy gap excitations polarized along the a and b axes, with gap energies 
A x = 0.48(1) meV and A 2 = 0.65(1) meV, respectively. These two gap values are 
consistent with the gap values we reported in Section 4.3. The experimental resolution 
is not sufficient to resolve these two modes. The gap in the c-axis polarized Haldane 
excitation, A 3 m 1.86 meV [77], is outside the energy range covered in the present 
experiment. In Fig. 4.13(b), at H = 5.5 T, one clearly sees the splitting of the 
doublet. Excitations polarized along the b axis (upper branch) are S z = states, 
z being chosen along the applied field. To first order, they are not affected by the 
magnetic field. The gap for the a-axis polarized Haldane excitation (lower mode) 
decreases with field by virtue of its mixing with the c-axis mode, induced by the 
Zeeman effect. 

The data were analyzed using a parameterized model cross section function, nu- 
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Figure 4.14: Field dependence of the gap energies (a) and excitation intensities (b) 
measured in NDMAP at T = 2.5 K (symbols), (c) A constant-Q scan collected at 
H = 5.5 T (symbols). Solid lines and shaded areas are as described in the text. 
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merically convoluted with the calculated spectrometer resolution. The cross section 
was written in the single mode approximation described by Eq. 4.1 and Eq. 4.2 in 
Section 4.3.3. [77]. For each of the two excitation branches visible in our experiment 
the dispersion relation was written as (hco q ) 2 = A 2 + v 2 sin 2 (qc) (model 1). The gaps 
and intensities for each branch were adjusted to best fit the experimental data. At 
each field the entire data set with energy transfers between 0.2 meV and 1.6 meV was 
analyzed with a single set of parameters in a global fit. The spin wave velocity v was 
fixed at v = 6.9 meV, as determined from the H = data set. The model was found 
to describe the data measured at up to 7 T very well, with a residual x 2 < 1-5 in all 
cases. The obtained field dependence of the gap energies is plotted in Fig. 4.14(a). 
The corresponding dispersion relations are shown in solid lines in Fig. 4.13(a)-(b). 
Fig. 4.14(c) shows a constant-Q scan extracted from the 5.5 T data set (symbols) 
plotted over the profile simulated using the refined parameter values (black line). 
The shaded areas show partial contributions of each of the two gap modes. Simi- 
larly, Fig. 4.15(a)-(c) shows a series of extracted constant-^ scans for H = 5.5 T, 
plotted on top of the corresponding simulated profiles (red lines). The measured field 
dependence of gap energies follows the general trend [solid lines in Fig. 4.14(a)] ex- 
pected from a perturbation treatment of the Zeeman term [83]. In agreement with 
the perturbation analysis [blue line in Fig. 4.14(b)], the measured intensity of the 
upper mode is field independent [Fig. 4.14(b), cyan symbols] while the intensity of 
the lower mode decrease with H as it acquires polarization along the c-direction to 
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which our experiment is insensitive [Fig. 4.14(c), green line]. In summary, at H < H c 
the behavior of NDMAP is well described by first order perturbation theory and is 
similar to that previously found in the extensively studied Haldane gap compound 
Ni(C 2 H 8 N2)2N0 2 (C104) (NENP) [72]. 

We now turn to the data visualized in Fig. 4.13(c), collected at H = 9 T> H c . The 
gap in the upper mode has increased substantially, to about 0.93 meV, as compared 
to 0.63 meV at H = 5.5 T. Interestingly, below this gap, a weak but well-defined area 
of scattering extends all the way down to the lowest accessible energy transfer. This 
is not conventional critical scattering associated with the imminent phase transition 
as the experiment was done away from (h, 0.5, 0.5) where elastic scattering develops 
for T < T N . 

This first experimental observation of quasi-one-dimensional quasi-elastic neutron 
scattering for H > H c implies that individual spin chains in NDMAP are in a critical 
phase from which quasi-long-range order is likely to emerge on cooling even in the 
absence of inter-chain coupling. Closer examination of the data reveals important 
changes in the high field excitations that go beyond closing the gap in the neutron 
scattering spectrum. Applying the low field analysis procedure to the 9 T data yields 
a significantly worse fit (x 2 = 2.0) than was achieved for H < H c . The problem is 
mainly in the wave vector and energy dependence of inelastic scattering at low ener- 
gies. The discrepancy is particularly pronounced in the shape of constant-E scans, 
where the data shows peaks considerably broader than simulations based on param- 




Figure 4.15: Symbols: constant-E scans collected in NDMAP at H = 5.5 T< H c 
(a-c) and H = 9 T> H c (d-f). Lines are as described in the text. Arrows indicate 
the projection of the resolution ellipse onto the chain axis. 
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eters obtained in the fit. To quantify this behavior we scaled the measured intensity 
by the energy transfer and integrated the result over the range 0.25-0.5 meV. The 
resulting ^-dependence was fit to a Gaussian profile to obtain an intrinsic (resolution 
corrected) width of 0.048(5)c*. The peak simulated using parameters determined in 
the global fit has an intrinsic width of only 0.023c*. The anomalous g-width is also 
visible in each measured constant-^ scan individually, as for example at 0.4 meV 
energy transfer [Fig. 4.15 (d)]. It is important to emphasize that the discrepancy is 
not related to an error in resolution calculation. Indeed, at lower fields, e. g. at 
H = 5.5 T [Fig. 4.15 (a)] Gaussian fits to const-E scans and simulated profiles yield 
identical widths. Moreover, even at H = 9 T at higher energy transfers, where the 
dominant contribution is from the 0.93 meV mode, the model cross section fits the 
data quite well [Fig. 4.15 (e) and (f), solid lines]. 

4.5.4 Discussion 

At least two theoretical frameworks could potentially account for the q-broadened 
quasi-elastic scattering in the high field phase of NDMAP. A classical easy plane 
antiferromagnet with a field in the easy plane can be mapped onto a Sine-Gordon 
model with a gapped spin wave spectrum and topological soliton excitations. [131] 
While solitons have a finite rest mass, they are stable topological objects that move 
freely along the spin chain. The soliton gas yields quasi-elastic neutron scattering 
with a q-width that varies in proportion to the soliton density. The soliton model 
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was used successfully to account for the phase diagram of NDMAP with fields normal 
to the chain axis. [119] Soliton pinning at defects may also account for the unusual 
quasi-two-dimensional frozen state in NDMAP for fields perpendicular to the easy 
plane. [76] Further evaluation of the Sine-Gordon model for the high field phase of 
NDMAP would require measurements of the H/T dependence of the q-width for 
comparison to the predicted soliton density. The free-fermion model presented in the 
introduction describes a different type of topological excitations. In this model at 
H > H c the lowest energy excitations with momentum transfers along the chains 
near q\\ = ir/c are single-particle or single-hole states. These have a dispersion rela- 
tion as shown by the red line in [Fig. 4.10 (b)]. To compare the 9 T data to the free 
fermion model, we constructed an empirical dispersion relation for the lower mode 
[model 2] with q F as a new parameter. The dispersion was assumed to be linear 
outside the region ir/c — q-p < q\\ < 7r/c + qp, and to follow half a period of a sine wave 
within these boundaries. The prefactor for the sinusoidal part was chosen to give the 
correct spin wave velocity at q\\ = ir/c ± qp. As illustrated in Fig. 4.10 (b), a given 
value of q F corresponds to a "negative gap" A = —2q F cv/n. With q F = 0.03(l)c* 
(A = —0.3(1) meV) model 2 yields a considerably better agreement with experi- 
ment (x 2 = 1-3) than model 1. The refined gap energy for the upper mode and 
the "negative gap" for the lower branch are indicated with triangles in Fig. 4.14 (a). 
Simulations based on resulting parameter values are shown in magenta lines Fig. 4.15 
(d)-(f), and the dispersion relations are shown in a solid line in Fig. 4.13(c). While 
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at high energies where the upper mode is dominant the two models are virtually in- 
distinguishable, below 0.6 meV transfer model 2 does a better job at reproducing the 
observed peak widths [Fig. 4.15 (d)]. The free fermion model makes a very robust 
prediction for the field dependence of q F . The Fermi vector q F determines the number 
of S z = 1-carrying particles present in the ground state at H > H c and is therefore 
directly proportional to the uniform magnetization: M = gji-Q 2q F /c*, where M is the 
magnetization per spin. For NDMAP the magnetization curve has been measured ex- 
perimentally [132] and at 9 T M « 0.1// B /Ni 2+ . This gives q F = 0.025c*, remarkably 
close to the experimental value. It is important to stress that the free-fermion model 
is a rather crude approximation. In fact, at H > H c the system is to be described 
of interacting fermions, i. e., a Luttinger liquid (for a recent discussion see 
Ref. [129]). As a result of these interactions, there is no well-defined Fermi surface or 
even single-particle poles in the dynamic susceptibility. Incommensurate correlations 
are obscured: the equal-time spin correlation function is a broad peak at the com- 
mensurate q\\ = 7r point. It is probably for this reason that static long-range ordering 
in NDMAP[76] and the related compound Ni(C5H 14 N 2 ) 2 N3(C104) (NDMAZ) [133] is 
commensurate. In the Luttinger liquid model the low energy excitations are a con- 
tinuum. Experimental evidence for this behavior was recently obtained in studies of 
NDMAP [119] and NENP [134]. While the resolution of our experiments on NDMAP 
is insufficient to distinguish between continuum and single-particle excitations, the 
free-fermion model appears to be a good starting point for the data analysis. It can 
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reproduce the measured data quite well, and yields a self- consistent estimate for qp. 

In summary, our experiments demonstrate q-broadened quasi-elastic scattering in 
the high-field phase of a S = 1 Haldane-gap compound. Various models based on 
topological excitations with low energy phase fluctuations may account for the data 
including the Sine-Gordon model has been shown to account for the phase diagram 
and the free fermion model that explains the observed q-width. Further theoretical 
and experimental work will be required to identify the correct description of the 
critical high field phase of the anisotropic spin-1 chain NDMAP. 
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Chapter 5 

Impurity in a spin-1 Heisenberg 
antiferromagnetic chain, Y2BaNi05 

5.1 Introduction 

Hole-doped quantum spin systems exhibit a variety of novel phenomena as a func- 
tion of dopant concentration. One class of such systems are the cuprates, which be- 
comes high temperature superconductors on doping. The dominant electronic and 
magnetic properties of cuprates are associated with the copper-oxide (CUO2) planes. 
In their undoped states, the compounds exhibit antiferromagnetic long range or- 
der below the Neel temperature T N . On introduction of a few percent of holes, 
the antiferromagnetic long range order is rapidly destroyed leaving behind a spin- 
disordered state. The hole-doped Cu02 plane is an example of a quantum spin liquid 
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in two-dimension. As its one dimensional analog, the S=l one dimensional Heisen- 
berg antiferromagnet, Y 2 BaNi0 5 , which can be doped with carriers, has attracted 
much attention [64, 88, 135, 136, 137, 138, 139, 140]. Here we present an inelastic 
neutron scattering study of the one dimensional spin liquid Y 2 BaNi05, with a focus 
on excitations in the bulk spin-1 chain and excitations associated with doping. 

5.2 Y 2 BaNi0 5 

The crystal structure of Y 2 BaNi05 is orthorhombic, space group Immm with 
lattice parameters o=3.7648 A, 6=5.7550 A, and c=11.324 A at T=10 K [141, 142]. 
Ni 2+ ions, which are at the center of compressed corner sharing oxygen octahedra, 
form chains along the a axis as shown in Fig. 5.1. Each of the S=l Ni 2+ ions is 
coupled to its neighbors through antiferromagnetic super-exchange via shared apical 
2 ~ ions. With Y 3+ ions positioned between the chains, these chains are magnetically 
isolated. Substituting Ca 2+ for Y 3+ introduces hole carriers in the compound, and 
they mainly reside in the — O — Ni — O — Ni— chains. No metal-insulator transition 
takes place, but the DC resistivity falls by several orders of magnitude [64]. Thus, 
Y 2 _ x Ca x BaNi0 5 is a one-dimensional analog of the cuprates, where off-Cu0 2 -plane 
chemical impurities donate holes to Cu0 2 planes. 
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Figure 5.1: A schematic view of the crystal structure of Y 2 BaNi0 5 . Ni0 6 octahedra 
stack along the a crystallographic axis, which is therefore the chain direction. The 
figure is taken from V. Sachan et al, Phys. Rev. B 49, 9658 (1994). 



5.3 Experimental procedure 
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Inelastic neutron-scattering measurements were carried out on the MAPS time- 
of-flight Fermi chopper spectrometer, at the ISIS pulsed spallation neutron source 
at Rutherford Appleton Laboratory in the United Kingdom. Single crystal samples 
of Y 2 _ x Ca x BaNi05 were grown by the Floating Zone Technique by Drs. T. Ito, K. 
Oka and H. Takagi from Japan. For every doping concentration, x=0, x=0.04, and 
£=0.1, two crystals of the same composition were coaligned to produce a sample 
with a mosaic of about 1° and mass of 10 g, 6.65 g and 4 g, respectively. Fig. 5.2 
shows a picture of a coaligned sample. The samples were mounted with the a and b 
axes horizontal. We index momentum transfer in the corresponding reciprocal lattice 
Q = ha* + kh* + Zc*, and we use q = Q a = 2nh to denote wave vector transfer along 
the spin chain. A closed-cycle refrigerator was used to cool the sample to T=6.7 K. 
Measurements were done with incident neutron energies of E =30, 60, or 90 meV. 
The sample was oriented with the a (chain) axis perpendicular to incident neutron 
wave vector kj. Chopper S was set at a frequency of 300 Hz. Each measurement was 
made with a total integrated proton current on target of around 3000 fj,A h. 

The MAPS spectrometer works as follows: A high intensity pulse of neutrons is 
produced by spallation every 20 mS. The neutrons are brought to useful energies by 
a water moderator and monochromated by a fast Fermi chopper. This produce a 
monochromatic beam of energy E , ranging from 15 to 2000 meV, with an intrinsic 
energy resolution of approximately 2 % full width at half maximum of the incident 
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energy E . A massive slow chopper running at 50 Hz reduces the fast-neutron back- 
ground. The monochromatic beam is scattered from the sample, and the scattered 
neutrons are measured in a 16 m 2 position sensitive detector array, with 160,000 pixel 
elements. The detectors are located at a distance of 6 m away from the sample po- 
sition. Low-angle detector banks cover scattering angles from 3° to 20°. High angle 
detector banks measure the scattering at larger momentum transfers and extend over 
scattering angles from 20° to 60°. From the time of detection of a scattered neutron 
in a particular detector, both the energy transfer hco and momentum transfer Q to 
the sample are calculated. Due to the high efficiency, high sensitivity and low back- 
ground of MAPS spectrometer, we were able to collect high quality data over a large 
(Q,w) range over the seven days beam time allocation assigned to the experiment. 

5.4 Experimental results and data analysis 

5.4.1 Y 2 BaNi0 5 

Magnetic one-dimensionality causes the parent compound of Y 2 BaNiOs to be a 
spin liquid. The material has a macroscopically coherent quantum ground state and 
lacks static order. The ground state is a valence bond solid in which the two spin-1/2 
degrees of freedom, that make up each spin-1 unit, are engaged in singlet formation 
with right and left nearest neighbors [53]. Excited states are triplets formed by 
broken valence bonds that propagate along the chain. Since a finite energy is needed 
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Figure 5.2: Two cylindrical Y 2 BaNiOs single crystals were coaligned to produce a 
sample of total mass 10 g and a mosaic of about 1°. The cylinder axis coincides with 
the chain axis. The length of the crystals is approximately 35 mm. 
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to create a triplet at rest, the triplet band is separated from the ground state by 
a Haldane gap [122]. (A more detailed discussion on the Haldane gap and on the 
ground state of a S=l chain can be found in Section 1.3.) Fig. 5.3 provides a view 
of the entire triplet band for Y 2 BaNi05 collected with initial neutron beam energy 
£^=90 meV. The figure shows the magnetic neutron scattering intensity multiplied 
by energy transfer, as a function of wave vector transfer q along the spin chains, 
and energy transfer, huj. (Incoherent elastic scattering and inelastic non-magnetic 
scattering has been subtracted from data as background.) At 8 meV, the onset of 
an intense mode of excitation is seen, which corresponds to the threshold energy 
(Haldane gap) for breaking a valence bond in the macroscopic singlet ground state. 
Magnetic neutron scattering are found only when hu and q lie on a well-defined 
dispersion relation, which is a direct evidence for a coherent triplet propagations 
along the chain. The overall magnetic band width is around 62 meV. Using the 
dispersion relation, hcu(q) = \/A 2 + v 2 sin 2 q, and the numerical results A % 0.41J, 
v ~ 2.49J[83, 108], we estimate that the intra-chain exchange constant, J, is around 
24 meV. The vanishing intensity for q=n2n suggests an isotropic singlet ground state. 
The steep dispersion along the chain direction is examined in constant-energy cuts 
above the gap energy along the wavevector transfer axis near the ID AFM zone center, 
q = 7r. The data are shown in Fig. 5.4. The two spin wave peaks are clearly resolved 
at high energies. The peaks were fitted to Gaussians whose positions are summarized 
in Fig. 5.5. We fitted the data to the dispersion relation hco(q) = \/ A 2 + v 2 sin 2 q to 
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extract accurate values of the Haldane gap and spin wave velocity. The fitting result 
are shown as a solid line in the figure. The fit gives A=9.2(2) meV, and u=60.5(4) 
meV. Gap values reported in Ref. [65] are 7.5, 8.6, and 9.6 meV for polarizations 
parallel to the a, b and c axis, respectively. Thus, the average of the directionally 
averaged gap value, A = (A a + A fe + A c )/3 = 8.6 meV, is close to our result (as we did 
not resolve the anisotropy splitting of of the Haldane gap). The obtained dispersion 
relation is also plotted in Fig. 5.3 in solid line. 

5.4.2 Y 2 _ a; Ca ;c BaNi05 

What happens to the quantum spin liquid when we add carriers by Ca doping? 
Fig. 5.6 (b) shows the entire triplet band in the £=0.04, Y 2 _ ;c Ca x BaNi05 sample. 
Fig. 5.7 (b) gives a close up look of the spectrum in energy transfer range < hco < 45 
meV. Fig. 5.7 (c) provides a view of the triplet band of x=0.1 sample in the same 
energy transfer range. Compared to the excitations in pure sample shown in Fig. 5.6 
(a) and Fig. 5.7 (a), one can conclude that excitations in Ca doped samples follow a 
similar dispersion relation as pure Y 2 BaNi05. To study how the Haldane gap energy 
and spin wave velocity change with doping, constant energy cuts were extracted from 
the data above the gap energy. Fig. 5.8 and Fig. 5.9 show the cuts for £=0.04 and 
x=0.1 samples, respectively. Solid lines are Gaussian fits to the peak. The positions 
of the peaks are plotted in Fig. 5.10; for comparison, the peak positions for the 
pure sample are also shown. The lines are least square fits to the dispersion relation 
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0.0 0.5 1.0 1.5 2.0 

q (it) 

Figure 5.3: Overview of magnetic fluctuations in Y 2 BaNi05 at T=6.7 K. The initial 
neutron beam energy was 90 meV. The chain axis was perpendicular to the incident 
beam direction. 
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Figure 5.4: Constant energy cuts at hco=8, 22, 36, 50 meV for pure Y 2 BaNiOs. The 
solid lines are described in the text. 
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Figure 5.5: Peak positions extracted from constant energy cuts at various energy 
transfers for Y 2 BaNi0 5 . The solid line is described in the text. 
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Figure 5.6: Overview of magnetic fluctuations in (a) pure (b) 4 % Ca doped Y 2 BaNi0 5 
at T=6.7 K. The initial neutron energy was 90 meV, and the chain axis was perpen- 
dicular to the incident beam direction. 
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Figure 5.7: The medium energy excitation spectrum of in (a) pure (b) 4 % and (c) 10 
% Ca doped Y 2 BaNi0 5 at T=6.7 K. The initial neutron beam energy was 60 meV, 
and the chain axis was perpendicular to the incident beam direction. 
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huj(q) = a/A 2 + v 2 sin 2 q. The fitting results are A=9.4(l) meV, w=62.2(l) meV 
for the £=0.04 sample, and A=9.9(2) meV. v=62.3(3) meV for the x=0.1 sample. 
Fig. 5.11 shows the extracted gap energy versus the hole concentration. For an 
isolated chain with L Ni 2+ spins, the energy of the lowest-lying one-magnon state is 
approximately [60] 



where A(oo) is the value of the gap for an infinite length chain. For low hole density, 
one would expect that the gap value will increase as A(L) ~ A(oo) + t4/L 2 . In £=0.04 
sample, the average chain length is L k, 1/x 25; In x=0.1 sample, the average 
chain length is L 10. Thus the gap value should increase with hole concentration 
as observed in the experiment. 

Fig. 5.12 shows the low energy detail of magnetic excitation in pure, £=0.04 
and £=0.1 sample. As opposed to the pure sample which clearly has no measurable 
intensity below the Haldane gap, there is considerable spectral weight below the gap in 
Ca doped samples. For the £=0.04 sample, Fig. 5.12 (b) clearly shows the appearance 
of two lobes of scattering displaced symmetrically about q = it. When we increase the 
concentration of holes to £=0.1 , these modes tend to spread out in energy as shown 
in Fig. 5.12 (c). Fig. 5.13 shows a cut through the sub-Haldane-gap feature along the 
wave vector axis at huj=2 meV and huj=5 meV. An "incommensurate" double-peak 
structure is clearly observed. Charge ordering could yield such an incommensurate 
spin order [143, 144, 145]. In the simplest picture, the hole density has the same 




(5.1) 
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Figure 5.9: Constant energy cuts at hco=8, 22, 36 meV for 10% doped Y 2 BaNi05. 
Lines are as described in the text. 
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Figure 5.10: Peak positions extracted from constant energy cuts at various energy 
transfers for pure, 4% and 10% doped Y 2 BaNi0 5 . The solid is as described as in the 
text. 
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Figure 5.11: Extracted gap values versus hole concentration x. 
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Figure 5.12: Low-energy details of the magnetic excitations in (a) pure (b) 4 % and 
(c) 10 % Ca doped Y 2 BaNi05 at T= 6.7 K. The initial neutron beam energy was 30 
meV, and the chain axis was perpendicular to the incident beam direction. 
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£=0.04 


x=0.1 




8q(n) 




8q(iv) 


k(h) 


2 meV 


0.033(1) 


0.08(3) 


0.045(2) 


0.09(1) 


5 meV 


0.041(1) 


0.07(2) 


0.046(2) 


0.10(3) 



Table 5.1: 8q and k values derived from fitting the data of Fig. 5.13 by the sum of 
two identical Lorentzians centered at 7r ± 8q with half widths at half maxima k. 

period as the spin density squared, so that 8q(x)/ir ~ x, where 8q(x) is the shift of 
the peaks from n. To explore the possibility, we fit our data to the sum of two identical 
Lorentzians centered at n±8q(x) with equal half widths at half maxima k. The results 
are listed in Table 5.1 and shown in Fig. 5.13 in blue lines. The results show that 
the shift 8q more or less remains the same while Tiui changes from 2 meV to 5 meV, 
suggesting the sub-gap mode is almost dispersionless in this energy range. Moreover, 
more than double doping concentration does not yield a proportional increase in the 
shift of peaks; rather 8q remains independent of x, which rules out the hole-ordering 
scenario. Another possibility is that Quasi-particle Quasi-hole pair excitations in 
Luttinger liquid produce an incommensurate magnetic fluctuations [135]; however, 
Sq(x) ~ x again is expected. Thus, our result contradicts this framework as well. 

Since the double peak structure is unchanged as we more than double the hole 
concentration, hence it should be interpreted as a single impurity effect rather than 
a result of correlations between impurities. We thus consider the magnetic structure 
factor for a single static hole donated to the NiO chain by the Ca ions. The hole 
is located on the superexchange mediating oxygen orbitals and carries a spin. It 
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q (7i) q fa) 

Figure 5.13: q cuts through the incommensurate peaks at Tiuj=2 meV and hu=5 meV 
for pure, 4% and 10% doped Y 2 BaNiOs. The dashed lines are the fits to the data of 
the pure sample. The solid lines are as described in the text. 
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induces an effective ferromagnetic (FM) interaction between the Ni spins on both 
sides [146, 147] shown in Fig. 5.14 (b). The FM interaction yields a ground state with 
inversion symmetry about the oxygen site. Furthermore, the spin density modulation 
developed around hole is not just confined to Ni spins neighboring the hole, but 
extends with an AFM modulation to a distance controlled by the correlation length 
of the quantum spin liquid, k -1 (See Fig. 5.14 (c)). The modulus squared Fourier 
transform of the spin polarization measured by neutron scattering for such a single 
FM bond inserted in an infinite chain can be written as [88] 

= wg f , (5 , 2) 

cosh k + cos q 

Eq. 5.2, which describes a function with pairs of peaks with half widths and incom- 
mensurability of order k, could explain the incommensurate double peaks. However it 
also predicts nodes at q = (2n + l)7r, which are not observed in the experiment. Thus, 
a more sophisticated model is required to fully account for our data. We first take into 
account the fact that we have a finite hole densities rather than an isolated impurity. 
Thus the polarization clouds developed around the hole overlap. A simple model 
is to truncate the polarization clouds at neighboring impurity sites. The structure 
factor, Si> n i(q), for a spin polarization cloud around n (n > 1) continuous impurities, 
truncated at /' spins to the left side, and / spins to the right side, involves the squared 
Fourier transform of the real space spin configuration and includes the contribution 
from the spin polarization cloud extending to the right and left side of impurity, from 
the nearest neighboring Ni spins to the impurity and the relative strength of the spin 
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of the hole on the oxygen site. The intensity we observed in the neutron scattering 
experiments is an overall average of the structure factors of all possible spin clusters 
around impurities, which can be written as 

S(q) = Y J Pi'niSvni{q)- (5.3) 

I'nl 

Here P\i n i = x n+2 (l — x) l+l ' +2 is the probability of finding such a spin polarization 
cloud, and x is the concentration of impurity. The summation in Eqn. 5.3 for S(q) 
can be carried out in exact analytical form: [87] 



S(q) = {2x 2 (l - x)\MM~ + ; "I!" + "'It. T ~"™ 9 ) 



'x 1 — e cosh k' + cos q 

x(l - x) _ _ 1__ 

1 — cos q v ' 1 1 - xe'i J 

q 2(1 -x) cosf + {cos q-x)S Q 

+S (S + 2 cos - - 1 - x) — 

2 1 — 2 cosg + x 2 

+ 2(l-x) 2 Re{Ml (1+ , , ) 2 . } 

+ 4s(l-s)i2e{-^(l + -^C^ T ) 
1 - xe % i 1 + e l, ?- K 

(^ + 5o e^ 2 )}}|F(Q)| 2 (5.4) 

where So is the relative strength of the spin of the hole on the oxygen site, \Moo\ = 
1/ (l+exp(iq — k)), and F(Q)| 2 is the magnetic form factor of Ni. Moreover, although 
holes are not fully mobile, they are delocalized over several lattice spacing as suggested 
by the transport measurement [64]. Eq. 5.4 is generalized to take into account hole 
spins distributed with exponential decay over / lattice sites. As shown by the black 
lines in Fig. 5.13, the model fits the data quite well. Assuming the hole spin has the 
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same net amplitude with either of the Ni 2+ spins next to the FM bond, the fit gives 
I ~ 5, k _1 =8.9(2) lattice spacing for £=0.04 sample and n~ l = 8.0(4) lattice spacing 
for #=0.1 sample. These values are comparable to the exponential decay length of 
6.03 lattice spacing calculated for the Heisenberg S=l chain. 

The data also shows that there is a fine structure in the excitation spectrum for 
Ca doped sample. This is most clearly appreciated in the cut along the energy axis 
at q = 7r shown in Fig. 5.15. In contrast to the clean gap in pure Y 2 BaNi05, which 
has no intensity below ~ 7 meV energy transfer, extra intensity is clearly observed 
at hco 5 meV for £=0.04 sample. An additional hump of intensity appears at 
Tiu) Ril2.5 meV. The origin of the resonance might relate to the composite spin degrees 
of freedom induced around holes. Finite interactions break the eight-fold degeneracy 
of the ground state associated with two chain ends spin-1/2 of Ni 2+ ions and a spin- 
1/2 of the intermediate hole. To figure out the energy values, we consider a simplified 
Hamiltonian 

H = JSi ■ S 2 + JS 2 • S 3 + J'Si • S 3 , (5.5) 

where J is the effective strength of coupling between the spin-1/2 of the chain end 
Ni 2+ ions and the spin 1/2 of the intermediate hole, and J' is the strength of the 
effective coupling between the two spin-1/2 Ni 2+ ions (\J'\ <C \J\). If the coupling 
between the spin 1/2 of the chain end Ni 2+ ions and the spin 1/2 hole is ferromagnetic, 
the lowest energy spin configuration has all the three 5=1/2 spins parallel. The 
eigenvalues of the Hamiltonian, ordered from low to high energies, are E = | + ^ 
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Figure 5.14: (a) Schematic of S=l/2 degrees of freedom at the edges of AKLT state, 
(b) Ferromagnetically coupled chain ends S=l/2 degrees of freedom (c)An Antiferro- 
magnetic droplet with a central phase shift of n and a size controlled by the correlation 
length of the quantum liquid. 
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(four-fold degenerate), E\ = IZ ^- (two-fold degenerate), and E 2 = —J + ^ (two-fold 
degenerate). The gap values are E x - E = -J' - |, and E 2 - E = The 
resonance at Tiu> ~ 5 meV is possibly due to the anisotropy splitting of the quartet, 
while the resonance at huj m 12.5 meV may arise from the the transition from the 
quartet to the lower energy doublet (Ei — E = — J 1 — |), which corresponds to 
J m 25 meV, taking into account that J 1 <C J. In the case of the antiferromagnetic 
coupling, the spin configuration is such that the two spin- 1/2 of Ni 2+ ions are parallel 
and the spin-1/2 of the hole is antiparallel to them. The eigenvalues are E = — J + ^ 
(two fold degenerate), E\ = — ^f- (two fold degenerate), and E 2 = |(2J + J') (four 
fold degenerate). Thus, the extra intensity at 5 meV could come from the anisotropy 
splitting of the lower energy doublet, and the intensity at 12.5 meV may correspond 
to the doublet to doublet transition {E\ — E = J — J'). Considering that J' is very 
small, we conclude that J ~ 12.5 meV. In £=0.1 sample, the features are smeared 
out in energy. In the sample with high hole concentration, rather than isolated spin 
polarons, they interact with each other via overlapping AFM polarization clouds, and 
produce effective coupling of random signs. Consequently, the peak in the energy scan 
is broadened. 

Another effect of doping is observed when we perform a constant q cut through the 
top of triplet band. Fig. 5.16 shows the energy cut at q=7r/2. The solid lines are the 
Lorentzian fits to the peaks. In pure sample, the peak is well defined. This indicates 
the energy cost to flip the spin which lies in the bulk of the chain is well defined. In the 
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Figure 5.15: Energy cuts at q = tt for (a) pure (b) 4% (c) 10 % doped Y 2 BaNi0 5 . 
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Figure 5.16: energy cuts through the top of the triplet band at q = 7r/2 for (a) pure 
(b) 4% doped Y 2 BaNi05. The solid lines are as described in the text. 

£=0.04 sample, however, the peak is broader, which reflects the the inhomogeneity 
experienced by a given spin on the chain due to the presence of impurities. It gives 
us a qualitative impression of the distribution of local environment seen by a single 
spin. The peak is broadened by about 1 meV in £=0.04 sample compared to the pure 
sample, which corresponds to a change in the exchange constant 8 J ' j J ~ 4.7%; this is 
consistent with the change in the spin wave velocity Sv/v = (62.2 — 60.5)/60.5 ~ 2.8% 
obtained from the analysis of the low energy part of the dispersion relation. 

5.5 Conclusion 

In summary, magnetic neutron scattering have been used to study the magnetic 
fluctuations in a doped one dimensional spin liquid Y 2 _a ; Ca ;! ;BaNi05. The parent 
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compound is a spin liquid with a cooperative single ground state and a Haldane gap, 
A=9.2 (2) meV. in the magnetic excitation spectrum. Replacing the off-chain Y 3+ 
by Ca 2+ introduces holes primarily onto apical oxygen atoms in the chains. We have 
shown that at energies above the Haldane gap the triplet excitations of the parent 
compound persist with doping. When we increase the hole concentration, an increase 
of the Haldane gap values is observed, with A=9.4(l) meV for £=0.04, and A=9.9(2) 
meV for x=0.1. Below the gap, calcium doping creates new excitations with charac- 
teristic wave vectors, which are displaced from the zone boundary by an amount of 
the order of the inverse correlation-length in the bulk spin chains. The incommensu- 
rability saturates with increased doping. The incommensurate peaks arise because of 
the characteristic spin density modulation that develops around an impurity in the 
singlet ground state of a quantum spin liquid. Our study on Y^zCa^BaNiOs gives a 
quantitative view of dilute spin polarons in a quantum magnet. 
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Chapter 6 
Conclusion 

In this dissertation, I have described an experimental study on three quantum 
magnets. Although they are all quantum spin chains, they exhibit a variety of fasci- 
nating cooperative phenomena. 

Specific heat measurements performed on the spin-1/2 chain CDC reveals a rich 
field-temperature phase diagram in this material. They also show an opening of a 
spin excitation gap in the magnetic field, which could be a result of an effective stag- 
gered field due to the alternating g-tensor and the Dzyaloshinskii-Moriya interaction. 
Such a spin-1/2 chain with a staggered field can be mapped onto a quantum sine- 
Gordon model, which is exactly solvable and from which many physically important 
quantities can be calculated. Elastic neutron scattering experiments have been per- 
formed to determine the magnetic structure of the ordered phase and the critical 
exponents of the ordering transition. Future inelastic neutron scattering experiments 
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can be performed to explore the rich excitation spectrum of the sine-Gordon model 
experimentally. 

NDMAP, a uniform spin-1 system, has an easy-plane anisotropy and forms chains 
along the hard axis. Zero field inelastic neutron experiments have been performed 
to characterize the spin Hamiltonian and determine the anisotropy split Haldane 
gap values. Elastic neutron scattering experiments show that NDMAP can only 
be magnetized in a field that exceeds the critical field at which the Haldane gap 
closes. Inelastic neutron scattering in a magnetic field below and above the critical 
field were performed to explore the field induced incommensurate soft modes which 
have been observed in the high field phase of a spin-1/2 chain. We found quasi-elastic 
scattering with the magnetic field applied in the easy plane and at temperatures above 
the phase transition to static order. However, the incommensurate spin correlations 
should occur in a gapless phase which is only expected to exist when a high field 
is applied perpendicular to the easy plane. Thus experiments need to be performed 
with magnetic field applied along the chain axis, while wave vector transfer has a 
significant component along the chain direction. Unfortunately, this geometry is not 
compatible with the conventional configuration where the field is perpendicular to the 
scattering plane. Thus advances in high field magnet systems for neutron scattering 
are needed to fully explore the high field phase of NDMAP. 

The spin-1 chain, Y 2 BaNi05, provides a unique opportunity to study the effects 
of doping by chemical substitution in a quantum spin liquid. Doping Y 2 BaNi0 5 
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with Ca introduce holes that reside on super-exchange mediating oxygen atoms. At 
energies above the spin gap, the triplet excitation of the parent compound persist 
with doping. Below the gap, we have observed interesting changes to the scattering 
cross section. Most importantly, magnetic scattering appears below the Haldane gap 
with an incommensurate structure factor. Our data reveal that an antiferromagnetic 
droplet with a central phase shift of n can account for the data. They also provide 
information on the internal energy levels of the droplet. Future neutron scattering 
experiments can be performed in magnetic fields to fully determine the level scheme 
of the spin polarons induced by bond impurities. 

Quantum spin chains, which display a broad spectrum of complex quantum me- 
chanical phenomena, as demonstrated in this thesis, provide a unique opportunity to 
test existing theoretical results and discover new cooperative behavior in many body 
systems. By accumulating knowledge of these relatively simple model systems, we 
can advance our understanding of the more complex problems, such as those involving 
mobile charge carriers. 
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